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Abstract 



Q_i The causal approach to perturbative quantum field theory is presented in detail, which goes back to 

^ a seminal work by Henri Epstein and Vladimir Jurko Glaser in 1973. Causal perturbation theory is a 

I . mathematically rigorous approach to renormaUzation theory, which makes it possible to put the theoretical 

setup of perturbative quantum field theory on a sound mathematical basis. Epstein and Glaser solved this 
^ problem for a special class of distributions, the time-ordered products, that fulfill a causality condition, 

C\| which itself is a basic requirement in axiomatic quantum field theory. In their original work, Epstein and 

^ Glaser studied only theories involving scalar particles. In this review, the extension of the method to theories 

^ with higher spin, including gravity, is presented. Furthermore, specific examples are presented in order to 

highlight the technical differences between the causal method and other regularization methods, hke, e.g. 
O dimensional regularization. 
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1 Introduction 

Quantum field theory (QFT) is more singular than quantum mechanics. The basic mathematical objects of 
quantum mechanics are square integrable functions, whereas the corresponding central objects in QFT are 
generalized functions or distributions. A potential drawback of the theory of distributions for physics is the 
fact that it is a purely linear theory, in the sense that the product of two distributions cannot consistently be 
defined in general, as has been proved by Laurent Schwartz [1], who was awarded the Fields medal for his 
work on distributions in 1950. 

If one is careless about this point the well-known ultraviolet (UV) divergences appear in perturbative 
quantum field theory (pQFT). The occurrence of these divergences is sometimes ascribed in a qualitative 
manner to problematic contributions of virtual particles with "very high energy", or, equivalently, to physical 
phenomena at very short spacetime distances, and put forward as an argument that the quantized version of 
extended objects like strings which are less singular than point-like particles should be used instead in QFT. In 
view of the fact that UV divergences can be circumvented by a proper treatment of distributions in pQFT, this 
argument for string theories is no longer convincing. 

We illustrate the problem mentioned above by a naive example of a "UV divergence" by considering 
Heaviside-6- and Dirac-5-distributions in 1 -dimensional "configuration space". The product of these two 
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distributions Q{x)S{x) is obviously ill-defined, however, the distributional Fourier transforms 



V2^J^{d}{k) = V2Tid{k) = J dx5{x)e-'^'' = 1, (1) 



^Q{k) = lim / dxe{x)e-'''''~''' = lim 

e\0 J e\0 



e\o k — ie 



k-iO' 



(2) 



exist and one may attempt to calculate the ill-defined product in "momentum space", which formally goes over 
into a convolution 

T{e5}{k) = ^ / dxe-"'''Q{x)5{x) = / dxe-"''' [ 4^e(A;')e+^'='^ / ^8{k")e+'^"\ (3) 
y2'K J yl-n J J y2'K J \/27r 

Throughout this paper, we use the symmetric definition of the (inverse) Fourier transform according to Eq. (7) 
and Eq. (8). Since / dx e*(*'+'="-'=)=^ = 2T:6{k' + k" - k), we obtain 

The obvious problem in x-space leads to a "logarithmic UV divergence" in k-space. It will become clear 
below that a concise description of the scaling properties of distributions, related to the wide-spread notion 
of the superficial degree of divergence of Feynman integrals, is crucial for the correct treatment of singular 
products of distributions in pQFT. 

In pQFT, the role of the Heaviside 9-distribution is taken over by the time-ordering operator. The well- 
known textbook expression for the perturbative scattering matrix given by 



n=0 

— oo — oo 



= XI ^ ? / "^^^1 • • • / d'^XnT[Hint{xi) . . . Hint{Xn)], (5) 
n=0 

where the interaction Hamiltonian Hint (t) is given by the interaction Hamiltonian density Hint (x) via Hint (t) — 
J d^x Hint{x), is problematic in the UV regime (and in the infrared regime, when massless fields are involved). 
A time-ordered expression a la 

T[Hint{.xi) . . . Hint{.Xn)\ = ^ ^{.x^n^ - ^nj . . . 0(a;n(„_i) - x^^jHint{.xn^) . . . HintixuJ (6) 

Perm. 11 

is formal (i.e., ill-defined), since the operator-valued distribution products of the Hint are simply too singular 
to be multiplied by 6-distributions. 

In this review, the construction of pQFT is reviewed from a causal point of view with a special focus on the 
regularization of distributions. Typical examples are discussed in the causal framework and compared to the 
corresponding treatment in the Pauli-Villars regularization or dimensional regularization. In the last section we 
describe a modern approach to quantum gauge theories including gravity. This shows that the gauge principle 
in a suitable formulation is a universal principle of nature because it determines all interactions. Therefore, 
any regularization method must be in accordance with it. 
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2 Mathematical Preliminaries 



2.1 Regularization of Distributions 

Distributions are continuous linear functionals on certain function spaces. There exist different spaces of distri- 
butions. For quantum field theory the most important function space is the Schwartz space 5(R"). It consists 
of infinitely differentiable complex-valued functions of rapid decrease, that means the functions together with 
their derivatives fall off more quickly than the inverse of any polynomial. The reason for the importance of 
S (R") is the fact that the Fourier transform (the expression px denotes a generalized n— dimensional Euclidean 
or Minkowski scalar product depending on the respective situation) 

= hp) = (27r)-"/' j f{x)e-'P^rx, (7) 

is a linear bi-continuous bijection from 5(R") onto »S(R"). Indeed, the inverse Fourier transform is given by 

:F-'{g}{x) = g{x) = (27r)-"/^ J gipY^^d^. (8) 

The dual space of »S(R") denoted by 5'(R"), is the space of tempered distributions. A tempered distri- 
bution d{f) is a continuous linear complex-valued functional on 5(R") 3 /; we also write d{f) — {d, /). 
The functions / G iS(R") are called test functions. The Fourier transform of a tempered distribution d is now 
simply defined by its action on the test functions: 

(^{4,/)'='(d,^{/}), /g5(R"). (9) 

In this way, by operating on test functions, various operations on distributions like differentiation, convolution 
etc are defined. Note that the definition Eq. (9), d{f) = d{f), is sometimes written in an intuitive manner by 
the help of formal integrals 

J d{q)f{q) dy = J d{q)f{q) d^ , (10) 

E" E" 

showing the close relation of the definition above to the Plancherel theorem. Of course, the common physical 
distinction whether the integration variable q is in "real space" or "momentum space" is of no relevance here. 

For mathematical details concerning the properties of distributions, we refer to [2, 3]. 

The most important distributions for field theory are related to Unear partial differential equations, for 
example the Klein-Gordon equation 

n-l 



An important distributional solution of the 4-dimensional Klein-Gordon equation is the Jordan-Pauli distribu- 
tion 

Dmi^) = ^ I d'p5(p' - m')sgn{po)e-'P^ (12) 

where the integral must be understood as a distributional Fourier transform; the factor i makes Dm{x) real. 
If we decompose the sign-function, sgn pq = Q{po) — Q{—po), we obtain the decomposition of Dm{x) — 
Dm^ {x) + Dm^ [x) into positive and negative frequency parts, for example 

i^W(^) = / - m')e(-po)e'^^ = J^^I ^'^'^^^^ " m')e(po)e-''^ (13) 
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In addition to these solutions of the homogeneous Klein-Gordon equation we need weak solutions of the 
inhomogeneous equation 

(n + m2)D™(a;)=5(a;). (14) 
This retarded distribution which vanishes for negative time ct = xq < (c is the velocity of light) is given by 

D^Jx) = e{xo)D"^{x) (15) 

and the corresponding advanced distribution by -D^(a;) = D"^^{—x). Finally, the so-called Feynman propaga- 
tor is defined as 

(^) = DZ(x) - Di;;\x) = DZ(x) + i^W(x). (16) 



Its Fourier transform is equal to 



„-ipx 

D^{x) = -{2n)-' I d'p-^ — -, (17) 



where the symbol iO stands for ie and the limit e ^ in the distributional sense, i.e. in <S'(R"). 

In the appendix, a concise list of the fundamental free field commutators and propagators is given, where 
we also explicitly account for the most common conventions concerning the signs and normalizations of the 
distributions. 

In standard QFT the Feynman propagator D'^ is associated with the inner fines of a Feynman graph in the 
simplest case of scalar particles (spin 0). In a lowest-order loop graph there arises the problem of multiplying 
two Feynman propagators D'^^{x) ■ D'l''^{x), a product which is ill-defined. In fact, in momentum space this 
product corresponds to a formal convolution of the form 

S(p) =C J d'kDmMDm,{p-k), (18) 

where C is a numerical constant; we shall always use the symbol C for constants which we do not compute 
explicitly because they are not interesting for our purpose. To simplify the notation the various C's stand for 
different constants. By counting powers of we see that the integral Eq. (18) is logarithmically divergent in 
the ultraviolet regime \k\ — > oo. To make it well defined we use a regularization of the Feynman propagator 

Df{x) 

nreg^^^ ^ ( ^ p - 

where C denotes a real normafization constant which depends on specifically chosen conventions. Modifying 
the Feynman propagator according to Eq. (19) at a high mass or energy scale given M is the basic essence of 
the so-called Pauli-Villars regularization. Note that the propagator term containing M has the "wrong sign" 
and does not correspond to the contribution of a heavy physical particle. However, for M — > oo, Dp^{k) 
converges to Dp{k) in the sense of tempered distributions. We present here one possible approach to calculate 
the scalar self-energy diagram. Using the Fourier transform 



1 



oo 

1 J ^is(k^-mHiO)^g (20) 







^2 _ _|_ 

the regularized propagator is equal to 

oo 

D';^{k) = C f ds e''^'-''' (e-''""" - e"^^^') . (21) 
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Substituting the Feynman propagators in the self-energy integral Eq. (18) by regularized ones, we obtain a 
finite integral 



oo oo 



E'^^g(p) J d^k j dsi j ds2e''''''-''°(^{e-'''"'' - e-'''^'^ x (22) 







Here the 4-dimensional /c-integral can be carried out by means of the Gauss-Fresnel integral 

j ^i{ak-+bk)^4j^ ^ ^expf-^V a> 0. (24) 
J la^ \ AaJ 



The result is 

no oo 

e-(si+S2)0 /iSiS2Si + S2 



oo oo 

2 



J J {si + s2y V p /V / 



(25) 







x(^e-''2™2 _e-*^2M2j_ ^^^^ 

Now we introduce the new integration variables ti — si/ {si + S2) and ^2 = 51 + S2, then we have 

1 00 

J:f^{p) J dtij dt2^e-'''+'''''^^-''^P' (27) 
00^ 

We have written ie for iO and take the limit £ ^ later on. To perform the i2 -integration we need the integral 

00 

/ Je— +--'^'J„(^) (29) 

a 

where we the lower limit of integration is a > in order to avoid the singularity at x = 0; note that Eq. (28) is 
integrable at t2 = 0. By differentiating twice with respect to z the denominator x"^ is removed and the integral 
can be easily evaluated 

piza—ea 

J'M ^ — . (30) 
Now the limit a — > is possible and two integrations in z yield 

Jo{z) = -iz[log{iz - s) - 1]. (31) 
Using this result in Eq. (28) the regularized self-energy integral becomes 



1 

^r^(p) = C dti zi \og{izi - e) - Z2 log{iz2 - e) z^ log(i^3 - e) + Z4 log(i^4 - £) 



(32) 



where 

zi ^ ti{l - ti)p^ - timj - {1 - ti)ml (33) 
Z2 = ti{l-h)p^ -timl-{l-ti)M^ (34) 
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^3 = ii(l - ti^ - hM' - (1 - ti)ml (35) 
^4 = ti(l - h)p^ - tiM^ - (1 - ti)M\ (36) 

The integral Eq. (32) still diverges for M — oo. We have to split off the divergent part. This process, 
called renormalization, must always be combined with regularization. In order to obtain a unique finite result 
we proceed as follows. We compute the special value E^''s(0) and subtract it from (16). Then the limit 

lim lim (Er^(p) - E-s(O)) "^E'ip) (37) 

is finite. It satisfies the normalization condition 

E'(0) = 0. (38) 

We will not calculate the finite self-energy explicitly because later on we shall discuss methods which 
give the result in a more elegant way. The subtraction of a constant in Eq. (37) is equivalent to the subtraction 
of a local term ~ 6{x) in x-space. If we considered a more singular distribution, then a certain polynomial in 
p must be subtracted which corresponds to a sum of derivatives of the 5-distribution in x-space. 

2.2 Scaling Properties of Distributions 

From the last section it is clear that the singular behavior of a tempered distribution d{p) at infinity or of 
its (inverse) Fourier transform d{x) at a; = is of central importance in QFT. To study these properties the 
so-called quasi-asymptotics of a tempered distribution is very useful. The definition is the following: 

Definition 2.1. The distribution d{x) G iS'(R™) has a quasi-asymptotics do{x) at x = with respect to a 
positive continuous function p{6), 5 > 0, if the limit 

\imp{5)5"'d{5x)^do{x)^0 (39) 

5— »0 

exists in 5'(E"). 

In the smeared out form of Eq. (39) with a test function (p e <S'(R'") 

hmp{5)(^d{x), cpQ)^ {do, cp). (40) 
we go over to momentum space to find an equivalent condition for the Fourier transform d{p). Since 

{d{x), vp(^))= {d{p), {vi^)y{p))=s-{d{p), m)) = {dQ, m), (41) 

where (p denotes the inverse Fourier transform, we get the following equivalent definition: 
Definition 2.2. The distribution d{p) e S'{K!^) has quasi-asymptotics do{p) atp = oo if 

\\mp{5)(d{^^, p{p))= {do, Cp) (42) 

exists for all ^ G 5(7?™). 

In momentum space the quasi-asymptotics controls the ultraviolet behavior of the distribution. Let us 
consider a scaling transformation 

Yimp{8){d{l), (p{ap)) = (io(p), (p{ap)) (43) 
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Since 

hmp{ad){d{^), m) = (Mp), Hp)) (45) 
5^0 ad 

exists, we may conclude that the limit 

lim ^ = .-rnJMphmi ^^p^^a) (46) 

exists, too, assuming that the denominator is different from 0. By another scaling transformation it follows 

Po(a6) = Po(a)po(&), (47) 

which implies po(ci) = o!^ with some real u. We therefore call p{5) the power-counting function. 

With help of the power-counting function we can now define the singular order of a distribution. 

Definition 2.3. The distribution d G S'{R^) is called singular of order uj, if it has a quasi-asymptotics dQ{x) 
at X = 0, or its Fourier transform has quasi-asymptotics do{p) at p = oo, respectively, with power-counting 
function p{5) satisfying 

lim^ = a", (48) 

for each a > 0. 
Eq. (46) implies 

a™(clo(p), ^(ap)) = (rfo(-), Hp)) = a-^{Mp), Hp)) (49) 

(X 

= {do{x), ^{-)) = a^'idoiax), cp{x)) = a-'^{do{x), ip{x)), (50) 
d 

i.e. do is homogeneous of degree cu: 

cio(-) = a-'^do(p), (51) 
a 

doiax) = a-^'^+'^^oix). (52) 

This impUes that do has quasi-asymptotics p{6) — 5'^ and the singular order u, too. A positive measurable 
function p{5), satisfying Eq. (48), is called regularly varying at zero by mathematicians [4]. The power- 
counting function satisfies the following estimates: If £ > is an arbitrarily small number, then there exist 
constants C, C and 5q, such that 

Cr+^ > p{5) > C'S""-', (53) 

for 6 < 5o. 

We want to apply the definitions to the following examples: 

1) d — 1: From Eq. (39) we get p{6) = 5~'^ and u — —m. 

2) d{x) = D'^Six) where 

(J f gai+...+am 

' + ■■■ + am- (54) 

OXi . . . ox^ 

Since 

dip) = {27r)-"'/'{tp)\ (55) 

we obtain p(6) = ^l"' from Eq. (42) and u = \a\. 

3) Let us consider the Jordan-Pauli distribution Eq. (12) which has the following form in x-space 

.Or mQ{x'^^ 



Dm{x) = — b{x^) - -^^J^{myJx^) , (56) 
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where Ji is the Bessel function. We shall write the ^-distribution with argument always in contrast to the 
positive scaling factor 6. The one-dimensional ^-distribution satisfies 



S{6'x') = (57) 



5^ 

whereas the term with the Bessel function stays bounded for SVx^ 0. Hence 



lim S'D^iSx) = ^-^5{x') = D,{x) (58) 

which is just the mass zero Jordan-Pauli distribution. This illustrates the general fact that the quasi-asymptotics 
do is given by the corresponding mass zero distribution. Since the Jordan-Pauli distribution is typically con- 
sidered in R"^ (m = 4), we find p{5) = and oj{Dm) — —2. 
4) The positive frequency part Eq. (13) 

D^r:^\p) = ^^{P'W - m") (59) 
is best considered in momentum space. Since 

j Q(^^y(^^-m^^ip{p)d*p = 6^ j e(po)5(p'-(5W)(^(p)d^p (60) 
^S^ [ , "^'^ + (61) 



we find 



limr^DW(0=I)W(p) (62) 



which implies oj{Dm^) — —2, in agreement with the foregoing example. We obviously have uj{Dm^) — —2, 
too. 

We notice from example 2 that the degree of singularity at x = increases with a; > 0. The distributions 
with negative u; have only mild singularities. This difference will be important in the next section. 



2.3 Splitting of Distributions 

In QFT the problem arises of multiplying certain distributions which are singular at x = by the discontinuous 
step function ©(xq), where Xq is time. We will consider this problem only for distributions with a causal 
support: Let d{x) — d{xi, . . . Xn) £ S'{BJ^) where xj e R^, m — An, be a tempered distribution depending 
on n space-time arguments. By 

F+(0) ^ {x\x'^ ^ xl - x'^ > , xo > 0} (63) 
we denote the closed forward cone, and by 

F^(0) = {x I > , xo < } (64) 
the closed backward cone. The n-dimensional generalizations are 

(0) = {(^1, ■■■Xn)\ Xj e F^(0) , Vj = l, . . . n}. (65) 
The distribution d{x) has causal support if 

suppdcr+(x)ur„(o). (66) 
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This means that all n space-time points are either in the forward light-cone or in the backward cone. The 
splitting problem now is to decompose such a distribution into a retarded minus advanced part 



d{x) — r{x) — a{x), (67) 

where supp r C r+ and supp a CT^. 

The simplest example of a causal distribution is the Jordan-Pauli distribution given by Eq. (56) where the 
splitting into retarded minus advanced distributions is trivially possible, see Eq. (15). It is misleading that the 
Feynman propagator Wl'ix) is also called "causal" sometimes, because it does not have a causal support due 
to the presence in Eq. (16). In the general case we have to distinguish two cases: 

a) Singular order a; < 0: In this case, the power-counting function goes to infinity by Eq. (53) 

p{5) ^ oo for 5^Q. (68) 

This implies 

We choose a monotonous C°° -function xo over with 

'O for t<0 
Xoit) = {<l for < i < 1 (70) 
1 fox t> 1. 

In addition we choose a vector v — (I'l, . . . Vn-i) € T"*", which means that all four- vectors Vj are inside the 
forward cone V^. Then 

n-l 

V ■ X — Vj ■ Xj — Q (71) 

i=i 

is a space-like hyperplane that separates the causal support: All products Vj ■ Xj are either > for a; e r+ or 
< for a; e r~. Then as a consequence of Eq. (69) the limit 

limxof — —]d(x)'^Q(v-x)d(x) — r(x) (72) 

<5^0 \ / 

exists. This is the case of trivial splitting where the multiplication by step function is possible. The result is 
independent of v. 

b) a; > 0: Now the power-counting function satisfies 

^^oo for 5^0. (73) 
To get a vanishing scaling limit as in Eq. (69) we choose a multi-index b with |6| = uj + 1 and consider 

{d{x)x\ ^(-)) = (ci(5|/)/,^(|/))5™+-+i ^ (doiy), ^ 0. (74) 

It follows that the splitting as in case a) is possible if the test function cp satisfies 

D^ipiO) = for \a\ < uj. (75) 
To achieve that, we introduce an auxiliary function w{x) e <S(R"*) with 

w(0) = 1 , D"w;(0) =0 for 1 < \a\ < cu, (76) 
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and define 

\a\=0 ' \b\=ui+l 

The function w{x) serves for the purpose of getting rapid decrease for \x\ —>■ oo. Now the decomposition 
according to a) Eq. (72) is possible 

{rix),if)''^{e{vx)d,Wip), (78) 

a{x) = r-d. (79) 

After construction r{x) defines a tempered distribution with suppr C r+(0). It agrees with d(x) on r+(0) \ 
{0} in the sense of distributions, because a test function e S with supp if C r+(0) \ {0} vanishes at a; = 0, 
together with all its derivatives, so that the additional subtracted terms in Eq. (77) are 0. But without these 
terms, there is no splitting of d{x) which makes sense for arbitrary (p e S, because the limit Eq. (72) exists 
on subtracted test functions only. If one does the splitting incorrectly by simple multiplication with Q{v ■ x) as 
in a), one is punished by the well-known ultraviolet divergences in quantum field theory. As we will discuss 
in detail later, these divergences appear in loop graphs which have a; > 0. For those graphs the naive splitting 
with Q{v ■ x) is impossible and, as a consequence, the Feynman rules do not hold. 
Again we have 

u{r) ^u{d) ^u{a), (80) 
This is a direct consequence of the definitions Eq. (77) and Eq. (78), because the limit 

limp(5)(r(x), ^{^))= ^mp{5)(d{x), Qwi^^(^))) (81) 

= \imp{5){d{x), {&W^)(^)) = {do{x), {QW^){x)) (82) 

exists with the same power counting function as d{x). But in sharp contrast to case a), the splitting b) is not 
unique. If f{x) is the retarded part of another decomposition, then the difference 

r-r=J2 CaD''6{x) (83) 

|a|=0 

is again a distribution with point support. Since c<j > 0, this time the splitting is only determined up to a finite 
sum of local terms according to Eq. (83). These undetermined local terms are not fixed by causality, additional 
physical normalization conditions are necessary to fix them. 

Before we proceed, it might help to provide some intuitive understanding of the distribution splitting pro- 
cess. One should remember the fact that the distributions appearing in local quantum field theory are more 
singular than ordinary functions, such that the products of the distributions are not necessarily well defined ab 
initio. E.g., the Feymnan propagator can be calculated in configuration space [5] 

^f(^) = - -^e{x')[MmV^^) - iN^{m^/^)] + ["^ Q{-x')K,{mV^), (84) 

47r 8nVx^ (^-kYsJ -x^ 

where J\,N\, and K\ are Bessel functions. For x^ ~ 0, Eq. (84) can be decomposed according to 

D'^ix) = + Oim'x') = Dlix) + 0(mV). (85) 

A formal product like D]^{x)Dp(x) contains the highly singular (formal) expression — iO)^ ~ l/x'^, 

and it is not trivial to understand the precise meaning of such a singular expression in the vicinity of the 
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light-cone where ~ 0. In general, the definition of distributional products works better in momentum 
space, where the analytic behavior of distributions appearing in quantum field theory is smoother and where 
the product goes over into a convolution. Still, the true difficulty is located in the point a; = 0, where the 
distributional behavior of the product of two Feynman propagators is no longer mathematically meaningful. 
In momentum space, this problem leads to a logarithmically divergent integral. This is the point where the 
causal method provides the well-defined tools to isolate this ill-defined part of the product from the regular 
part on R''\{0}. Generally, perturbation theory itself is unable to describe local, "zero-distance" interactions 
without further input. At least, if the mathematics is done correctly and distributions are treated correctly, then 
all results remain finite at every calculational step, i.e. well-defined. The point x = is essential because 
the distributions are most singular there. The subtracted terms have no direct physical meaning because they 
remain with free parameters, this is the freedom of (finite) renormalization, which is discussed in further detail 
in the sequel. 

For practical reasons, explicit calculations in quantum field theory are usually done in momentum space. 
As a natural consequence, we must investigate the splitting procedure in p-space. We need the distributional 
Fourier transforms 

T-'[Q{v ■ x)]''^ m (86) 
T-^[x''w\{p)^{iD^Yw{p). (87) 

Since 

(D»(0) = {{-TD% ip) = {-r{D^S, ip) (88) 

= (-)»(27r)-/^((-ip)", ^) = (27r)-/^((ip)", (89) 
we conclude from Eq. (78) that 



|a|=0 



(90) 



= (27r)— /^(x*^^>-E-- - /' (91) 

|a|=0 

where the asterisk means convolution. We stress the fact that the convolution x * d is only defined on subtracted 
test functions, not on alone. Interchanging p' and p in the subtraction terms, we may write 

(f, (^) = (27r)— dkx{k){d{p-k)-{2n)-^/^J2^-^P'' J dp' d{p' - k)D;,w{p'), <p) ^. (92) 

a 

After partial integration in the p' -integral this is equivalent to the following result for the retarded distribution 

a r 

dk m [d{p -k)- (27r)— /2 J] ^ y {D^Ap' - ■ (93) 

|a|=0 

Here the A; -integral is understood in the sense of distributions as in (59). 

By considering the Fourier transform of Eq. (83) we see that f(p) is only determined up to a polynomial 
in p of degree u. Consequently the general result for the retarded distribution reads 

u> 

rip) = rip) + CaP" (94) 
l«l=o 
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with f{p) given by Eq. (93). We now assume that there exists a point q e R"* where the derivatives D^f{q) 
exist in the usual sense of functions for all 1 61 < lo. Let us define 



|6|=0 



(95) 



This is another retarded distribution because we have only added a polynomial in p of degree uj. Furthermore, 
this solution of the splitting problem is uniquely specified by the normalization condition 



D^rq{q) = 0, l&l < u. 



We compute 



D^f{q) = (27r)-"/' / dkxik) liD^iq - k) 



(96) 



(97) 



a—h 



b<a 



a — h)\a\ 



dp' w{p')D^,d{p' - k) 



from Eq. (93) and substitute this into Eq. (95). Since 



E 

b<a 



{p-qf q"-^ 



b\ {a-b)\ a\ 



b<a 



the subtracted terms in Eq. (93) drop out 



f,{p) = (27r) 



-m/2 



dkx{k) dip — k) 



\b\=0 



(98) 



(99) 



(100) 



This is the splitting solution with normalization point q. It is uniquely specified by Eq. (96), that means it does 
not depend on the time-like vector v in Eq. (86). The subtracted terms are the beginning of the Taylor series 
atp = q. This is an ultraviolet "regularization" in the usual terminology. It should be stressed, however, that 
here this is a consequence of the causal distribution splitting and not an ad hoc recipe. 

It is well-known that causality can be expressed in momentum space by dispersion relations. Therefore 
we look for a connection of the result Eq. (100) with dispersion relations. We take q = in Eq. (100), 
which is possible if all fields are massive, for example, and consider time-like p E T^. We choose a special 
coordinate system such that p = (p?, 0, 0, . . .). Note that this coordinate system is not obtained by a Lorentz 
transformation from the original one, but by an orthogonal transformation in R"*. Furthermore we take v 
parallel to p, i.e. v — (1, 0, 0, . . .). Then v varies with p, but this is admissible because Eq. (100) is actually 
independent of v. We now have Q{v • x) — ©(x?) and the Fourier transform (54) is given by 



X{k)-i'^^r^'~'S{kl,k2,...kn-l] 



k^ + iO' 



We always use the mathematical notation iO for ie with the subsequent distributional limit e 
result in Eq. (100) we shall obtain 



(101) 
0. Using this 



ro(p?) 



2^ 



k^ + iO 



d(p°-^?,0,...) 



(102) 
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-E^(-mV(??-^?,o,... 

a=0 ' gO=0- 



(103) 

The transformation of this result to the usual form of a dispersion integral leads to the following result: 



+00 

—00 

The proof is given in [6], proposition 2.4.1. This expression is a subtracted dispersion relation. To write down 
the result for arbitrary p e r+, we use the variable of integration t — kg /pi and arrive at 

+00 

fo(p) = — / dt. (105) 



—00 



For later reference we call this the central splitting solution, because it is normalized at the origin (q — 
in Eq. (96)). The latter fact has two important consequences, (i) The central splitting solution does not 
introduce a new mass scale into the theory. If g 7^ 0, then \q'^\ = defines such a scale, (ii) Most symmetry 
properties of d{p) are preserved under central splitting, as we will see later, because the origin g = is a very 
symmetrical point. In the self-energy computation in sect. 2. 1 by means of Pauli- Villars regularization we 
have also calculated this central solution Eq. (38). 

It is easy to verify that the dispersion integral Eq. (68) is convergent for |i| — > 00. But it would be 
ultraviolet divergent, if lj in Eq. (105) is chosen too small. Consequently, the correct distribution splitting with 
the right singular order u is terribly important. Incorrect distribution splitting leads to ultraviolet divergences. 
This is the origin of the ultraviolet problem in quantum field theory. 

3 Perturbative S-Matrix Theory 

In perturbation theory all quantities are expanded in terms of free fields. To decide which free fields are 
relevant we notice that all interactions in nature can be described by quantum gauge theories, gravity included. 
Therefore, it is important to discuss quantized free gauge fields and their gauge structure. The latter is defined 
by means of ghost fields. In contrast to the functional integral approach to QFT where the Faddeev- Popov 
ghosts play indeed a somewhat ghost-like role, these are genuine dynamical fields in the causal approach. 
Regarding regularization it is a subtle problem to perform it in a way such that gauge invariance is conserved. 
In this respect dimensional regularization is technically advantageous. 

3.1 Free Fields 
3.1.1 Scalar Fields 

First we consider a neutral scalar field with mass m which is a solution of the Klein-Gordon equation 

(n + m2)(^(x) = 0. (106) 
A real classical solution of this equation is given by 

d^p 



^{x) = (27r)-=^/2 j ^ (a{p)e-'^^ + a*(p-)e^^^) , (107) 



where 



px = p^x^ — p- X = Pi^x'^, E = +\/jP + m? = p". (108) 
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In quantum field theory a{p) and a*{p) become operator- valued distributions satisfying the commutation rela- 
tions 

[aip},a\q)] = 6^''\p-q}. (109) 

where 5*^^) denotes the Dirac's 5-distribution, all other commutators vanish. The quantized Bose field is now 
given by 

ip{x) = (27r)-3/2y" ^ (^a(p)e-'^'" + at(p)e^*'"). (110) 

The cross denotes the Hermitian conjugate. In fact, there exists a Fock-Hilbert space representation of the 
Bose field which proves the consistency of the quantization, (p is obviously Hermitian (/5^(a;) = (/^(x). 

Let us call the second term in Eq. (1 10) involving the creation part 99*^+^ and the first term with a{p) the 
absorption part ip'-~\ Then by Eq. (109) their commutator is equal to 

[^^-\x),^^^\y)] = (27r)-^/2| ^e-^^(^-^) = (Ill) 
In the same way we get 

[ip^^\x),^^-\y)] = -iD^-){x-y). (112) 
Then the commutation relation for the total scalar field is given by the Jordan-Pauli distribution 

[ip{x), ip{y)\ = -iD^{x - y). (113) 

The charged scalar field involves a slight generalization of the neutral one: 

ip{x) = {2n)-''^^ J ^L(^a{p)e-'P'' + b\p)e'P''y (114) 

It contains two different kinds of particles whose absorption and emission operators satisfy 

[a{p),a\q)] = 6{p-q) = m,b\q)] (115) 
and all other commutators vanish. Then it follows 

[ip{x),ip{yy] = -iDrn{x-y) (116) 

but [ip{x),^{y)] = 0. 
3.1.2 Spin-1/2 Fields 

Spin- 1/2 fields are needed to describe leptons and quarks. Spinor fields are solution of the Dirac equation 

i-f^d^'il){x)^m'^{x). (117) 
The 7-matrices obey the anticommutation relation 

+ = Sy"''- (118) 
To define the quantized Dirac field we consider a solution of Eq. (1 17) of the following form 

^{x) = (27r)-=^/2 j d'p J2 Mp)e-'^''bM + Vs{p)e'^''dl{p)i (119) 
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The u- and v-spinors herein are obtained from the Fourier transformed equations 

(Pi^r - m)us{p) = (121) 

(P/.7^ + m)^;,(-p) = 0, (122) 

with the normalization 

u\{p)Us'{'p) = 6ss' = vl{p)v,i{p) (123) 

u\{p)Vs'{-p) = = vl{-p)Us'{p) (124) 

ul{p)-f°Us'{p) = ^Sss' = -Vs{p)^-f°Vs'{p). (125) 

The u- and I'-spinors span the positive and negative spectral subspaces of the Dirac operator, respectively, 
which are defined by the projection operators 

s 

P-iP) = T.<-P)''l(-P) = (^)^°' ^127) 

s 

where 

^ = PmT'', Po^E^ + m2. (128) 

The projections are orthogonal 

P±{p)' = P±{P), P+{P) + P-ip) = 1. (129) 

The quantization of the Dirac field is easily achieved by considering the b's and d's as operator-valued 
distributions satisfying the anticommutation relations 

{bs{p),bl,{q)} = S,.,S^'\p- q) = K(p),4(g-)}, (130) 

and all other anticommutators vanish. We do not treat the Majorana case for neutral spin- 1/2 fermions here. 
Then, the operators 6 and d can be interpreted as annihilation and their adjoints 6^ d'^ as creation operators 
and the Fock space can be constructed from a unique vacuum in the usual way ([7], sect. 2.2). To get the 
anticommutation relations for the whole Dirac field we need the adjoint Dirac field 

ij\x) = (27r)-^/^ j d'p [bl(p)Us(p}^e'P^ + dsip}Vs{p}'' e-'^% (131) 

Multiplying by 7°, we get the so-called Dirac adjoint 

^(x) = il}\x)f = ^^^^ + (132) 

= (27r)-=^/^ j d'p h\{p)uMe'''^ (133) 

i^^~\x) = {2n)-^/''J d'pd,{p)MP)e-''^- (134) 
With the aid of Eq. (130) we find 

{^l^i-\x), i^i'\y)} = (27r)-^ j d'pu,,(p)u,,{^e-'P^^-yl (135) 
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In the result Eq. (135), the covariant positive spectral projection operator Eq. (126) appears 



{^(-)(a;), ^^+^1/)} = (27:)-' j + m)e-^P^^-y^ -^S^+\x - y) (136) 

= ~i{ip + m)Dl+\x-y). (137) 
In the same way, one obtains the other non-vanishing anticommutator 



(138) 



^-i{ip + m)D^-\x-y). (139) 
This gives the anticommutation relation for the total Dirac field 

{'il^{x),i^{y)}^-Sm{x-y), (140) 

with 

Sm{x) = S^^\x) + = {i^ + m)Dm{x - y). (141) 

The anticommutators between two ip's and two ip's vanish. 

3.1.3 Vector Fields 

Next we consider the massless vector field which obeys the wave equation UA^^{x) = 0. Its Fourier decompo- 
sition reads 

A^{x) = (27r)-3/' j 4i fa^(^)e-'*^^ + a'^(^)te*'=") , (142) 
and it is quantized in Lorentz-invariant form according to 

\A^{x\ A^{y)\ = g^^iD^{x - y). (143) 
We need also the commutators of the absorption and emission parts alone 

[A^_{x), A\{y)] = g^''iD^^\x - y), (144) 

A''_{y)] = g^''tDi-\x - y). (145) 

We are working in the Feyman gauge for the sake of convenience and covariance. However, asymptotic 
massless spin 1 particles only have two polarization degrees of freedom. Consequently, the four polarization 
types of emission operators introduced above create unphysical particle states, and the space of physical states 
is a subspace of the full Fock-Hilbert space. 

This observation is closely related to the issue of gauge transformations, therefore we also comment here on 
gauge transformations in perturbative quantum field theory. The massless vector fields describe non-interacting 
photons and gluons in the standard model. In classical electrodynamics the vector potential can be changed by 
a gauge transformation 

A'l'ix) = A^ix) + Xd^uix), (146) 

where u{x) is assumed to fulfill the wave equation Du{x) — because we want the transformed field A'^{x) 
to satisfy the wave equation also. In QFT, the quantized A'^{x) should fulfill the same commutation relations 
Eq. (143) as A^{x). This is true if the gauge transformation Eq. (146) is of the following form 

A'^'{x) = e-'^^A^'{x)e'^^, (147) 
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where Q is some operator in the Fock-Hilbert space. Expanding this by means of the Lie series 



^A^'{x)-iX[Q,A''{x)]+0{X''). (148) 
and comparing with Eq. (146) we conclude 

[Q,A^{x)]^id^u{x). (149) 

The operator Q will be called gauge charge because it is the infinitesimal generator of the gauge transfor- 
mation defined by Eq. (146). For the following it is important to have Q nilpotent 

Q2 = 0. (150) 



The important consequence of this property is the fact that the factor space Fph = Ker Q/Ran Q is isomorphic 
to the subspace of physical states. Here, Ran is the range and Ker the kernel of the operator Q. The overline 
denotes the closure; note that Ran Q is not closed because is in the essential spectrum of Q. We will not 
discuss this in detail but refer to the literature [6]. Such a nilpotency according to Eq. (150) is characteristic 
for Fermi operators. Therefore, we assume u{x) to be a fermionic scalar field with mass zero, a so-called ghost 
field. This field has the following Fourier decomposition 

uix) = (27r)-=^/2 y ^ (^C2{p)e-'P- + Ci(p)te^p-) . (151) 
In addition, we introduce a second scalar field 



u[x) 



(27r)-3/2y" A^(^_ci{p)e-'P^ + C2{p)^e'P^y (152) 



The absorption and emission operators Cj, c\ obey the anticommutation relations 

{c,(p), Ck{q)^} = 5jj,5^^\p- q). (153) 

Some remarks are in order here. Firstly, also bosonic fields would do the job in the case of an abelian theory like 
quantum electrodynamics (QED). However, non-abelian gauge theories like quantum chromodynamics (QCD) 
require fermionic ghosts. In order to avoid any conflict with the spin-statistics theorem, states containing ghosts 
necessarily do not belong to the physical sector of the Fock-Hilbert space of the theory under consideration. 
Secondly, when the spin-1 fields become massive, as it is the case in the standard model for the W^- and Z- 
boson fields, the corresponding ghost fields also become massive. The ghost mass then depends on the chosen 
formalism (i.e., the gauge fixing, [8]). In the following, we allow the ghost fields to be massive, but m — 
holds true whenever the related vector fields are massless. The formaUsm used below is chosen so that the 
ghost mass is equal to the vector boson mass. 

Again, the absorption and emission parts (with the adjoint operators) are denoted by (-) and (-I-). They 
satisfy the following anticommutation relations 

{u^-\x),u^+\y)} = {2n)-'J ^ e^^^^-^) = -iDM{x-y) (154) 

{u^+\x)M'\y)} = -(2^)"'/ ^e^^(^-^) = -tDi^\x-y). (155) 
All other anticommutators vanish. This implies 

{u{x),u{y)} = -iD^{x - y) (156) 
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and {u{x),u{y)} — 0. Then it is not hard to verify that the nilpotent gauge charge Q satisfying Eq. (149) is 
given by 

g = ^ d^x [d^A'dQU - [pQd^A'y] = J £x d^A'^dou (157) 

where the integrals are taken over any plane x^ — const. 

Now we return to the defining property of Q as being the infinitesimal generator of gauge transformations 
given by Eq. (147) and Eq. (149). We introduce the notation 

dQF^[Q,F], 

if F contains only Bose fields and an even number of ghost fields, and 

dQF^{Q,f}^QF + FQ, (158) 

if F contain an odd number of ghost fields. Then dg has all properties of an anti-derivation, in particular the 
identity 

{AB,C}^A{B,C}-[A,C]B (159) 

implies the product rule 

dQ(F(x)G(y)) = (dQF(x))G(y) + (-1)^- F(x)dQG(y), (160) 

where nj? is the ghost number of F, i.e. the number of u's minus the number of u's. The gauge variations dg 
of our free fields now are 

dgA^" = id^u, dgA'^ = id^Ui (161) 

dgu = 0, dgu = {Q,u} = -id^A^", dgu± ^ -id^A'^. (162) 

The latter follows from the anticommutation relation Eq. (156). dg changes the ghost number by one, i.e. a 
Bose field goes over into a Fermi field and vice verse. Then the nilpotency Q'^ — implies for a Bose field Fb 

dgFB = {g, [g, Fb]} = Q{QFb - FbQ) + {QFb - FbQ)Q = 0, 
and for a Fermi field F 

dgF = [Q, {Q, F}] = Q{QFb + FbQ) - {QFb - FbQ)Q = 0, 

hence 

4 = (1.4.31) 

is also nilpotent. The gauge variation dg has some similarity with the Becchi-Rouet-Stora-Tyutin (BRST) 
transformation [33, 34] in the functional approach. However, the BRST transformation operates on interacting 
fields (mainly classical) and the quantum gauge invariance which we are going to define is completely different 
from BRST invariance. 

Now we consider massive vector fields. These fields will be used to represent the W^- and Z-bosons of 
the electroweak theory, for example. They obey the Klein-Gordon equation 

{U + m^)A^'{x) = Q. (163) 

Since a spin-1 field has three physical degrees of freedom, we need one subsidiary condition to define unphys- 
ical states. As this we can choose the Lorentz condition 

d^Ai'^x) = 0. (164) 
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The commutation relations are similar to the massless case, for example (see Eq. (143)) 



[A^^{x), A'^iy)] = g^^'tD^ix - y). (165) 

where only the massless D-distributions must be substituted by massive ones. 

As in the massless case we would like to characterize the physical subspace with help of a nilpotent gauge 
charge Q. The old definition given by Eq. (157) does not work because nilpotency is violated: 

. ^ ^j^^ ioufa^u ^ -i,^J^.u3^ ^ 0. (166) 

To restore it we modify the expression for Q by introducing a scalar field $(a;) with the same mass m as the 
gauge field A'^{x) 

Q = J d^x {d^A"" + m$) ioM. (167) 
All fields satisfy the Klein-Gordon equation 

(n + m^)$ = (168) 

(□ + m^)u = 0, (169) 
but, while u{x) is a Fermi field, is quantized with commutation relations 

[^{x),^{y)] = -iD^{x-y), (170) 

and all other commutators are the same as before. Now we can check the nilpotency: 

= -i j (fx [d^A" + m$, Q] = 0, (171) 

because the first term in the commutator gives —iDu and the second one —iw?u, so that the sum is zero by 
Eq. (169). The infinitesimal gauge transformations or gauge variations in the massive case are now given by 

dQA''{x) = [Q, A^'ix)] = id^'uix) (172) 

dg^ix) = [Q, = imu{x) (173) 
dQu{x)^{Q,u{x)}^0 (174) 

dQu{x) = {Q, u{x)} = -i{df,A^ + m$(a;)). (175) 

The last equation follows from Eq. (167); using Q"^ = from Eq. (175) implies Eq. (173). 

Let us stress the difference between our approach to massive gauge fields and the conventional one. In 
the usual approach one starts with massless gauge fields and the scalar field $ is the so-called Goldstone 
boson. The fields become massive after "spontaneous breaking" of the gauge symmetry. We start directly with 
massive vector fields. To define a gauge variation with a nilpotent Q, we are forced to introduce the scalar 
field $, spontaneous symmetry breaking plays no immediate role. There is a common misconception in the 
literature, that the Higgs field "gives mass" to the particles. One could also argue that if particles are massive, 
then a consistent theory requires additional degrees of freedom, i.e. a Higgs sector. 
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3.1.4 Spin-2 Fields 

Finally, we comment on spin-2 quantum gauge theories which can be analyzed on the same footing as spin-1 
theories. We only consider the massless case which is relevant for quantum gravity. We start from a symmetric 
tensor field h°'^{x) with arbitrary trace which is assumed to satisfy the wave equation 

= (176) 

The gauge transformation similar to Eq. (146) is of the form 

lap ^ ^aP ^ ^ ^/3,a _ ), (177) 

where the comma denotes partial derivatives. This transformation leaves the so-called Hilbert condition 
= unchanged, if fulfills the wave equation 

□m" = 0. (178) 

The Hilbert gauge condition is analogous to the Lorentz condition in the spin-1 case. The corresponding gauge 
charge can immediately be written down in analogy to Eq. (157): 

j d^xh^>',^doUa. (179) 

The vector field Ua must be quantized with anticommutators, in order to get Q nilpotent. The operator Q given 
by Eq. (179) is the right infinitesimal generator for Eq. (177) if it has the following commutator 

[Q, = (u"^f^ + M^'" - g^'^u^',^ ) {x) = - ib'^^^'^u^,^ = dgh""^. (180) 

The factor — i/2 is convention and the 6-tensor 

^a^M^ = K^"^^'^'^ + ^""^'^^ - g^'^g^l (181) 

often appears in connection with tensor fields. It is also the four-dimensional extension of DeWitt's supermetric 
[9]. The commutator Eq. (180) implies the following commutation relation for the tensor field 

[/i"^(a;), h^^'iy)] = -ib'''^^'' Do{x - y) (182) 

= -^(^"V^ + rg^^ - g^'^gnM^ - y)- (183) 

The vector field must again be quantized with anticommutators 

{u^{x),u'{y)} = ig^''D,{x-y) (184) 

and the anticommutators between two n's or two tt's vanish, u^^ are called vector-ghost fields, respectively. 
An explicit representations of the ghost fields is given by 



(185) 



u^{x) = (27r)-=^/^ ( + - g^^ di{p)^ e''^) , 

where the absorption and creation operators satisfy the commutation relations 

{c^{p),cl{ky} = 6,,6^^6^'\p- k). (186) 
The gauge variation of the vector ghost fields now follows from Eq. (179) 

dgu" ^{Q,u''}^0 (187) 
c?Q«" = {<5,it"} = i/i"^^. (188) 
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3.2 The Causal Structure of the Perturbative S-Matrix 
3.2.1 General Construction of the S-Matrix 

Perturbation theory reUes strongly on the axiom of causality, as shown by H.Epstein and V.Glaser [10] after 
previous work of Stuckelberg, Bogoliubov and Shirkov [5]. The S-matrix is constructed inductively order by 
order as a formal power series of operator valued distributions 

OO OO ^ „ 

S{g) = -^^"H^) = 1 + E 1 / ^'^1 ■ ■ ■ d^^nTnix,, . . . xn)g{xi) . . . g{xn) (189) 

n=0 n=l ^' 

where g{x) is a tempered test function that switches the interaction. The operator- valued distributions T„ act 
in the Fock space of some collection of free fields. They are called time-ordered or chronological products and 
should verify the so-called Bogoliubov axioms. 

(1) It is clear from Eq. ( 1 89) that T„ can be assumed to be completely symmetrical in all variables xi,...Xn. 

(2) We must have Poincare invariance: 

X\i ■ ■ ■ Xji )U-l = T^{A-xi + a,...A-x,, + a) (190) 

for all proper Lorentz transformations A. In particular, translation invariance is essential in the causal approach. 

(3) The central axiom is the requirement of causality which can be written compactly as follows. If X = 

{xi, . . . Xm} € K.''™ and y = {yi, . . . ?/„} G R''" are such that Xi — yj ^ for all i and j, we say X is later 
than Y, X > Y. We use the compact notation T„(X) = Tn{xi, . . . x„) and by X U F we mean the union of 
the elements of X and Y. In particular, the expression T„_,_^(X U Y) makes sense because of the symmetry 
property (1). Now the causaUty axiom expresses causal factorization: 

r„+„(xur) = r„(x)r„(F), vx>r. (i9i) 

Physically this means that later action does not influence what has happened before. 

Like S{g) given by Eq. (189), the inverse S{g)''^ can be expressed by a perturbation series 

S{g)-' = 1 + 5] ^ / ^^^1 ■ ■ ■ d^^nfnixi . . . xn)g{xi) . . . g{xn) (192) 

n=l 

The corresponding n-point distributions T„, called anti-chronological products follow from Eq. 
inversion of a power series 

n 

r=l Pr 

where the second sum runs over all partitions Pr of X into r disjoint subsets 

X ^XiU...UXr, Xj^iJ}, \Xj l^rij. (194) 

All products of distributions in Eq. (193) are well-defined, because the arguments are disjoint sets of points 
such that the products are direct products of distributions. 

(4) Unitarity of the S-matrix 

S{g)-' = S{gy. (195) 
can now be expressed by means of the time-ordered products in the form 

f,(X) = T„(X)t. (196) 



(189) as formal 
(193) 
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It is one aim of QFT to prove unitarity for the physically interesting theories. In the inductive construction 
unitarity is not used. 

Now we are ready to turn to the inductive construction of T„ (xi . . . a;„) starting from Ti (x) which is a given 
interaction Lagrangian or coupling. Suppose all Tm{xi,. . . , Xm) for 1 < m < n — 1 are known and have the 
above properties (l)-(3). Then, according to Eq. (193), the Tm{X) can be calculated for all 1 < m = |X| < 
n — 1. From this it is possible to form the following distributions 



A'^{x,...Xn)^J2^n,{X)Tn-n,{y,^n) (197) 

(X), (198) 

where the sums run over all partitions 

P2: {xi,...,Xn-i} = XUY, X^0 (199) 
into disjoint subsets with \X\ — rii > 1, \Y\ < n — 2. We also introduce 

...Xn)^ R'ni^l^ ■■■Xn)- ^n(^l' " " -^n)- (200) 

If the sums are extended over all partitions P2, including the empty set X = 0, then we get the distributions 

)^J2'^n, iX)Tn-n, (Y, Xn) (201) 

^A'^ + Tnixi...Xn), (202) 

Rn{xi, . . . Xn) = ^ Tn-m (Y, Xn)fn^ {X) (203) 

^R', + Tn{xi...Xn). (204) 

These two distributions A^, Rn are not known by the induction assumption because they contain the unknown 
Tn{xi, . . . Xji). Only the difference 

D^^R'^-A'^^R^- A„ (205) 

is known according to Eq. (200). What remains to be done is to determine Rn (or An) in Eq. (205) separately. 
This is achieved by investigating the support properties of the various distributions. 

We recall the definition Eq. (65) of the n-dimensional generalizations of the forward and backward light- 
cones with respect to the point x. Then we have 

supp Rn,+i {Y, x) C r+ +1 {x) (206) 

and 

suppAn,+,{Y, x) C r- +i(x). (207) 

Because of these support properties, R and A are called retarded and advanced distributions, respectively. The 
distribution D, which can be expressed by Eq. (205), then has a causal support: 

SUppDn{xi, . . .Xn-l,Xn) C 1+ (x^) U (Xn) • (208) 

We do not present the proof here (see [7], sect. 3.1) but we can indicate the essential reason for this important 
causal support property: According to Eq. (197) and Eq. (198) D„ is a sum of commutators 

D'^{X^ ...Xn) = Y,[Tn-n,{y,Xn), fn,{X)]. (209) 
P2 
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Since all T's are products of free fields, the commutators contain Jordan-Pauli distributions which have causal 
support according to Eq. (56). 

Now we see the inductive construction clearly before us: From the known Tm{xi, . . . , Xm), m < n — 1 
one computes An given by Eq. (197) and Rn from Eq. (198), and then Dn — R'^ — A'^- One decomposes Dn 
with respect to the supports Eq. (208) 

Xl ) • • • ) Xn ) = Rn{xi, ...,Xn)- An{xi, ...,Xn), (210) 

suppi?„ C r+_i(a;„) , suppA„cr~_^{xn). (211) 
Finally, r„ is found from Eq. (202) or Eq. (204) 

r„(a;i, ...,Xn)= Rn{Xu ...,Xn)- R'ni^U • • • , ^^n) (212) 

= An{xi, . . . - A'^{xi, . . .,Xn). (213) 

The only non-trivial step in this construction is the distribution splitting Eq. (210). In sect. 2.3 we have 
discussed the splitting of causal numerical distributions. The transformation of the operator- valued distribution 
Dn to numerical distributions is achieved by means of Wick expansion. Dn can be written in a unique way in 
terms of normally ordered products of free field operators ^ 

-Dn = 5^4i...fe„(a;i, . . -a^n) -^kiiXl) ...^kA^n) (214) 

where /ci, . . . , A;„ are indices specifying the individual field operator types and the corresponding field compo- 
nents related to external symmetries (e.g., Lorentz indices) and inner symmetries (e.g., color indices), depend- 
ing on the theory under study. 

In the normal product between double dots all absorption operators stand to the right of all emission 
operators. Consequently, the vacuum expectation value of a normal product vanishes. This allows us to 
write the normal ordering of, e.g., n scalar fields (p in the compact form 



'x^)...^{xn)= J2 (^,<^'"'H^i)---<^'"'"(^n)^^):<^'n^i)---<^'i^n):, (215) 



and the generalization to products of general free field operators is straightforward. Here O, is the vacuum in 
Fock space and the brackets (•, •) mean the scalar product. Note that all terms with an odd number of sj — 
are zero because the vacuum expectation value vanishes. 

The splitting of Dn in Eq. (214) can now be carried out by splitting all numerical distributions dk. We have 
learned in sect. 2.3 that the splitting solutions may be not unique. Then the free but finite local terms must be 
chosen such that all properties required for the S-matrix are true. This is a subtle problem for gauge theories. 



3.2.2 Example 

As a simple illustration of the causal method we consider the coupling 

Ti{x) = -iX : ip\x)(p{x) : ^{x) (216) 

between a charged scalar field (p of mass m and a neutral scalar $ of mass M. This theory has the same 
structure as quantum electrodynamics: (/? is a spin-0 electron and $ a scalar massive photon. To perform the 
normal ordering according to Eq. (215) we need the commutators, or contractions 

<R^\y) = [<^(x)(-),(^t(^)(+)] ^ (Q^ cp{x)cp^{y)) = -iDl+\x - y) {111) 

and similarly for $. 
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In the inductive step from Ti to T2 we must first compute 

^^2(^1,^2) = -ri(x2)ri(xi) (218) 

by normal ordering: 



= X^(-. (p\xi)(p{xi)(p\x2)(p{x2) :: ^Xi)^X2) : 

-iD^)(x2 - xi) : ip\xi)(p{x2) :: $(xi)$(x2) : 

-iDl^\x2 - xi) : <^{xi)^\x2) :: <^{xi)<^{x2) : 
-iD'-j^\x2 - Xi) : ip^{xi)(p{xi)ip'^{x2)(p{x2) : 
-^^(^^ - ^i)^L+Ha;2 - xi) : *(xi)*(x2) : 
-L'W(X2 - Xi)L'i]:)(x2 - xi) : ip\xi)ip{x2) : 
-L'^^(a;2 - Xi)d''^\x2 - Xi) : (y9(a;i)v?"^(a;2) : 

-zDW(x2 - - Xi)Di^\x2 - Xi)). (219) 

We emphasize that here the product of two D^"'")(x)D*^+)(x) is well defined, in contrast to the product of two 
Feynman propagators Df{x) in sect. 2.1. The reason is that in the Fourier-transformed expression of, e.g., 

j^2 J d\ D^r^^ iP - Q)D^r^^ (?) (220) 

the intersection of the supports of the two IDm'^ is a compact set. This can be easily understood if one re- 
members that the support in momentum space of the two individual Pauli- Jordan distributions in Eq. (220) is 
contained in a forward and a backward light-cone, respectively. 

In the same way A'2 can be computed and then D2 — R2 — A2. From the third line in Eq. (219) we get the 
following contribution to D2: 

Dl = -i\^Dm{x2 - xr) : ^{xi)lp\x2) :: ^{x^)^{x2) : . (221) 

(+) 



Here can be trivially spUt. The retarded part with respect to X2 contains D'^{x2 — x-i). Adding R2 ~ D 
and using D'^ + D^r^^ = we finally obtain 

Tl = tX'D^{x2 - X,) : ip{x,)ip\x2) :: $(xi)$(a;2) : . (222) 

This gives "electron-photon" scattering in this model. The one contraction between the vertices Xi and X2 is 
represented by the Feynman propagator. So in tree graphs the usual Feynman rules hold. 

Now let us consider a loop graph with two contractions, for example "vacuum polarization" which comes 
from the fifth term in Eq. (219). The corresponding causal distribution is given by 

Dj{x^,X2) = XWiy)Di;l\y) - D^^\-y)D^^\-y)] : ^{x,)^{x2) :, (223) 

where y — x-y — X2. We calculate 

d+{-y)'^ D^^\y)D^^\y) {llA) 

in momentum space: 

d+{k) = -(27r)-^y" d"pe(A;o -po)<^((A: -p)' -m2)e(po)<5(p' -m^). (225) 
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It is easy to evaluate this for time-like k in the special Lorentz frame such that k = (/cq, 0). The result for 
arbitrary k is then 



(226) 



The total result for the Fourier transform of the square bracket in Eq. (223) denoted by d{k) is simply obtained 
by substituting ©(/cq) by sga{ko). 

Since d{k) has a constant quasi-asymptotics it has singular order a; = 0. Consequently, the distribution 
splitting is non-trivial. Therefore, the two internal lines in the vacuum-energy graph cannot be represented 
by Feynman propagators; the Feynman rules are not true for loop graphs. Instead we must use the dispersion 
integral Eq. (104) with uj — 0: 



+ CX) 



f{ko) 



-kn 



<hh) 



27r J (po - «0) {ko - po + «0) 



(227) 



(27r: 



1-3 



' 1 - ^^sgn(po) 



(228) 



For the time-ordered product we have to calculate £(/co) —fiko) — f'{ko). r'{ko) ~ d+{—ko) nicely combines 
with Eq. (228) so that 

+00 

1 



t{ko) = 



tvn 



ds- 



(27r)4 4 J s(ko -s + iO) 

Am? 



Am? 



(229) 



This integral is elementary, the final result for time-like momentum k^ — \J—k'^ is equal to 

1 + ( ^ ~ \/ -'- ~ 4^ 1 

1 + 



t{k^) = 



(2vr) 



1 ^ 1 ., 

4m2 j ^ flko 



1 (2m aA 4r^) 



Am?) 



(230) 



The result for space-like k is obtained by analytic continuation. We leave it as an exercise to the reader to 
compare this result to the alternative outcome given below at the end of sect. 4.2. 



4 Regularization Methods 

4.1 Basic remarks 

It would not be worthwhile to recapitulate the well-known details of the different regularization methods which 
are on the market. While the Pauli-Villars regularization can be considered as an ad hoc method to solve the 
apparent problem of infinities in perturbative quantum field theory, dimensional regularization is treated by 
many introductory texts like [11, 12], including the original works on the topic [13, 14]. 

Dimensional regularization regularizes Feynman diagrams by analytic continuation to 4 — e (complex) 
space-time dimensions and isolates infrared and ultraviolet divergences as poles in e. From a technical point 
of view, the main question is how to evaluate Feynman diagrams in n dimensions, i.e. it is necessary to know 
the usual Feynman rules of the theory, properties of Dirac matrices in n dimensions, and techniques Uke the 
Feynman parametrization for performing the momentum integrals in n dimensions. Difficulties may arise 
when one has to deal with topological quantities which exist only in integer dimensions. 

The same applies to the causal method, however, there is a big conceptual difference between the causal and 
dimensional regularization approach to perturbative quantum field theory. Whereas the causal approach is a 
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mathematically fully understood perturbative method, is is hard to give a precise meaning to the idea of physics 
in an arbitrary complex number of space-time dimensions. Still, the method has many interesting technical 
advantages, and although a proof of the physical equivalence of the causal and dimensional regularization 
method is lacking due to the technical complexity of the problem, one should not be too pessimistic about that 
issue. 

In the forthcoming section, we will illustrate to conceptual differences by presenting some specific exam- 
ples in the light of the different approaches. There, the diagrams will turn out to be finite in most cases, but 
this is not the central issue since infinities can always be removed in one or the other way. In the present 
short section, we present a direct comparison of the treatment of the scalar one-loop integral in the causal and 
dimensional regularization approach. 



4.2 Scalar One-Loop Diagram in n Dimensions 

The positive frequency part of the Pauli- Jordan distribution in n dimensions is given in the causal framework 
by 

D^^\x) = J rpS{p' - m')e{po)e-'P'. (231) 

Strictly speaking, the expression above is well-defined in integer dimensions. In order to obtain an analytic 
expression for the scalar one-loop diagram in arbitrary dimensions, we generalize Eq. (225) to 



m 



^ I d-pe{k'-p')5{{k-pf-m')e{p')5{p'-m'). (232) 



(27r)3-/' 

We exploit the last 0- and ^-distributions in Eq. (232) and evaluate the integral over space-like momenta 

/jji-i 
^^S{kl - 2k^E)Q{k'^ - E) (233) 
2E 



and obtain a radial integral, using E = \p\ = y/kl/A — rn^ 

^^"^ = r("(n-l)/2) ®^^° " ^^'^ j ^l^^l^^(2^°(^72 - E))Q{k' - E), (234) 

since the (n — 2) -dimensional surface of an (n — 1 ) -dimensional unit ball is given by (27r('^~^)/^)/r((n — 1)/2). 
Substituting | = EdE, can be written 



r((n-l)/2)^ ^ ... ,^vw 2A;0 



or 

?rn— 3 



dfik) = ,3 TVl^rMk' - ^m')e{k Y ^^ ■ (236) 



(27r)3«/2-2r((n-l)/2)^''^ ^^v-w 
This result can be compared directly to Eq. (225) for n = 4. 
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Now, the interesting point is that the (1 — t + iO)-term in the central splitting formula Eq. (105) generates 
the real part of the scalar loop amplitude, denoted here by i'^^^k), since 

^ = P— i7r5(l - t), (237) 

where the symbol P denotes the Cauchy principal value in the sense of distributions. We obtain for k in the 
forward light-cone 

lHe(i>)(A;)) ~ e{k' - 4.m') , (238) 

where we have omitted numerical factors. In fact, this result is valid for arbitrary momenta k. 
Now the scalar loop integral in n dimensions is given by the expression 



XW(p) = f — (239) 

J {k^-m^ + iO)[{k-pf-m^ + iO] 

Feynman parametrization 

1 

= J da[aA + {l- a)B]~^ (240) 



and a subsequent momentum translation k^ ^-^ k^ + ap^ leads to the formal integral 

1 

I<".W = /«/.„[.^-™^ + .(l-„y + «]-. (241, 



To evaluate this integral, we may use the relation 

d^'f' ■ " r(2-f) ^^^^^ 

= « TT 2 / , (242) 



which is divergent for n > 4. In the present case, we have 

a^^m^- a(l - a)p^. (243) 

The calculation of the finite integral in the case n = 3 is discussed in detail in sect. 5.3. We focus here on 
the case n — A. Note that the following manipulations are formal to some extent. Since 

and 

a"-4 = l + (n-4)loga + ..., (245) 
where 7 is the Euler-Mascheroni constant, we obtain in the limit n — > 4 

1 



2i7r f 

T^^\p) — >• ivr^ / dalog(m^ — a(l — a)p^ — iO) — i7r^7 — ivr^ log(7r). 

4 - n J 


A relevant finite part of the integral above is given by 

1 

^itliP) = -^^'J dalog[^ ^-^^ j 



(246) 
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= da ( , "if-;', J . (247) 

J \m"^ — q;(1 — a)p'^ J 



where an integration by parts was performed. For > 4m^, the — iO-term generates a real part of Ireg{p), 
which is obtained in a trivial manner from the partial fraction decomposition 

— ^ — ^ + — = -2 — "2! 2 — 7^ (248) 

with 



11/ 4m2 



ai,2 = -±-^/l-— , 0<ai,2<l, (249) 



leading to 



mt{4:l{p)) («i - «2)e(p2 - 4m') = Q{p' - - (250) 



in accordance with Eq. (238). For > 4m , one can also write 



This solution is normalized according to 



]im4f(p) = 0, (252) 



i.e. it corresponds to the central splitting solution in the causal approach when continued analytically to 
arbitrary p. 

We observe that the real part of the scalar loop diagram coincides both for the causal and the dimensional 
approach, and it is straightforward to show that this result holds in arbitrary dimensions. Furthermore, Eq. 
(236) provides a kind of "dimensional" generalization of the causal method. From the real part, the imaginary 
part of the amplitude is obtained from the dispersive splitting formula in the causal approach or by direct 
computation according to the rules of dimensional regularization. Up to finite renormalizations, the finite parts 
of the results also agree. 



5 Comparison of Regularization Methods to the Causal Approach: 
Specific Examples 

5.1 Axial Anomalies 

Axial or triangle anomalies are a subtle problem because their treatment by regularization of divergent Feyn- 
man integrals is unsatisfactory. On gets the impression that the anomalies are a consequence of the ultraviolet 
regularization. Then the question remains whether by some other method of calculating the divergent integral 
the anomaly might disappear. The causal method is free of such uncertainties as we are going to show. 
We consider QED with pseudovector and pseudoscalar couplings 

Ti{x) ^ icvjv{x)Aij,{x) +icAjA{x)Bu,{x) 

+icJ^{x)U{x). (253) 



29 




Figure 1: Axial vector triangular graph Figure 2: Pseudoscalar triangular graph 

Here 

j!^ = : ^^/^^ : , = : ^^f^^^^jj ; (254) 
are the vector and axial vector currents and 

j^ = i : V^7^V^ : (255) 

is a pseudoscalar, all being formed from a free massive Dirac field %jj{x) with mass m. The vector-, axial vector 
and pseudoscalar vertices defined by Eq. (253) will be abbreviated by V, A and 11 in the following. The fields 
Af^, and Il{x) play no essential role and are, therefore, assumed as classical external fields. 
From Eq. (254) and Eq. (255) we have the following divergence relations for the free currents 

d^Jv = 0, J A = 2m j^, (256) 

as a consequence of the Dirac equation. Our problem is whether similar divergence relations hold at higher 
orders, in particular for the two triangular graphs with vertices VVA— and VVH of Figs. (1) and (2) which 
contribute to the 3-point function T3. To compute the latter, we must first calculate 

Ds{xi,X2,X3) = T2(xi,X3)T/(x2) + T2{x2, X3)TI{xi) + Ti (x3)Tj(xi , X2) 

-T/(xi)r2(x2, X3) - r/(x2)T2(xi, X3) - T^{xi,X2)Ti{xs), (257) 

where we have used unitarity to express the T-distributions of the inverse S-matrix. Concerning the triangle 
graphs, the 2-point distributions which contribute, come from Compton scattering. For these distributions the 
usual divergence relations in Eqns. (256) still hold, so that 

%d''J'''^'{xuX2,X3) = 2m^C^^(a:i,a:2,X3). (258) 



dx^ 

Here dB.d.,, are the numerical 3-point distributions corresponding to the triangle graphs displayed in Figs. 
(1) and (2) without the external fields A, B, 11. The question is whether the same relation remains true after 
splitting for the retarded distributions. One therefore defines the anomaly by 



The t-distributions have the same anomaly because the r'-distributions are anomaly-free 



= -7n^r''B^"''{xi,X2,X3)-2m^r>:^^''{xi,X2,X3). (259) 

OXn 
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Since we work with massive Feraii fields, we can perform the splitting in momentum space by means of 
the central solution 



hoo 



— oo 

where p, q are assumed to be in the forward cone V~^, and similarly for r^. a; is the singular order of the d- 
distributions. The Fourier transformation is carried out in the difference variables yi—Xi — x^, y2 — X2 — x^, 
taking translation invariance into account. From Eq. (258) we then get 

i{p,s + q,,)dT'''{P. q) = 2m^<i'^^(p, q), (261) 

and the anomaly in Eq. (259) becomes 

£^i^W2^p^ g) = i(p^^ + q^^)f%'>'^''^ - 2m— f^^'^^ (262) 

Substituting Eq. (260) and the analogous equation for f^^ herein, and using Eq. (261), we arrive at the following 
formula for the anomaly 

00 

—00 

Hence, the anomaly is due to the fact that — ujb 

To evaluate Eq. (263) we only need the pseudoscalar d-distribution. From the first three terms in Eq. (257) 
we find 

+i755i+)(a;3 - X2)r'S^-\x2 - x^)r'SF{x^ - x^T 
+^75-5^(2:3 - ^2)7''^ -51+^^2 - xi)7^i>Si-)(xi - X3)] + tr [xi ^ X2, III ^ 112] . (264) 

Here 5"^^ denotes the anti-Feynman propagator which is obtained from the Feynman propagator in momentum 
space by changing +iO into —iO: 

SUp) = (27r)-^ , (265) 

It comes from the adjoint in Eq. (257). If one replaces Sm^ by Sm^ and vice versa without changing the 
arguments, one gets a'^^^'\ The difference r' — a' gives d^^^^. 

Expressing the spinor distributions by scalar ones, we see that the terms with three p contain 75 plus five 
other 7-matrices. Then the trace vanishes. In the non-vanishing terms one has at least one factor m instead of 
p. This lowers cu by one, so that we get cj^ = 0, instead of the power-counting estimate 1. But we will see 
below that the splitting with = or 1 gives the same result. If we replace 275 by 7^^7,5 in Eq. (265), we 
get the r'-distribution for the axial vector graph: r''^'^^^^^. Then the terms with three ^ contain 75 plus six 7^* 
matrices and the trace does not vanish. In this case we have the power-counting result = 1. 

After Fourier transformation 

f'^^^'ip, q) - (27r)-^ j r'^'^'iyi, y2)e'^^'+'^^''^ dyi dy2 (266) 



we shall obtain 
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+il,S^{-P + k)^^-Sl+\-p + k)r'S^-\k)\ + tr [p ^ g, „ |. (267) 
Here we have introduced P — p + q. Computing the trace we get 

X {[I-{P,p) + I+{q, -p) + /+(p,P)] + b " (268) 
where the Lorentz invariant integrals /± are given by 

I±{P. q)= J d^kQ{-k'')S{k^ -m^)x 

X e{k' - p^iik - pf - rT?\- (269) 

[k — qY — m^±i\) 

Owing to the two (5- and 9-functions, these integrals vanish if p is not in the region p^ > Am?, po < 0. But 
if p is in this region, one can use JC\_ invariance 

/±(Ap,Ag) = /±(p,g), VA e 4, (270) 

to choose Ap = (— a/p^^ 0). Then the integration is done as follows: first we integrate over /c°, using the first 5. 

„ — * — * 

In the spatial integration d'^k, we use polar coordinates with q as polar axis. Integration over | A;| kills the second 
5, while the integral over the azimuth cp gives trivially 27r. The remaining integration over cos '& — k ■ q/ {\k\\q\) 
is elementary. The result in an arbitrary Lorentz system is equal to 

I±{p,q) = ^e(V)e(p^- 4m2)-^x 



f-pq + g2 + ^(1 _ 4myp^)N ± tO\ 

X log , , (271) 

y-pq + q^- ^y{l-4:m^/p^)N±iQ^ 

where 

N = N{p, q) = {pqf - p'^q^. (272) 
The final result for Eq. (268) is now given by 

r'^'^'i.P, q) = s'^^^^'V^/^r (p, g), (273) 
r'(p, q) = JI^^^ [_e(_/)e(/ _ 4^2) ^^g^ 



where 



(27r)VA^ 

-e{-q')e{q' - Am') log, -e{-P'')e{P' - Am') log,] , (274) 

^ ^ rl+m±V^^^±^l (275) 

yq-'+pq- ^y{l-Amyp^)N + iO^ 

fp' +pq+^{l- 4myq^)N + iO \ 

log2 = log I \ o / o =Y^ ^ ] ' (276) 

\p2 j^pq- ^(1 _ Am^/q^)N + 
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( -pq+^(l-AmyP^)N-iQ ^ 

logo = log , . (277) 

V-pg - ^(1 - 4mVP2)iV - iQJ 

The expressions for a' and d are similar: 

^/MiM2 ^ et'^i'^^f^p^qpa'ip, q), (278) 



-e(g°)e(?2 - iog2 -e(p°)e(p2 - 4m2) logg] , (279) 

JmM2 ^ ^MiM^a/?^^^^ J^^^ ^) ^ (280) 



mc?/C. 



(^(P) = ""'"^ sgn(p")9(p^ — 4m^) log^ 



(27r)VA^ 

+sgn(g°)e(g2 - 4m^) log^ +sgn(P°)e(p2 - Am^) log,] . (281) 
Since the scaling limit in Eq. (281) is equal to 

]imd^^'^-(Xp,Xq)^C!t^(p,q), 



A— »oo 



we conclude that cUj, = 0. However, the central splitting solution is independent of choosing a; = 1 or a; = 0, 
respectively. To see this, we calculate the difference 



' 01=1 ' u)=0 



f,.,—n — 



oo 

27r J 1-t + iO t) 

— oo 



oo 



~27r J ~ ^' 

— oo 

because the denominator is an odd function of t. 

To calculate the anomaly, we now insert Eq. (280) into Eq. (263) 

d^'i'^p, q) = £'^i^^-'^p„?^a(p, q), (282) 

where 



oo 



XTTl f dt 

a{p,q)^—CyCA / —d{tp,tq), (283) 



t 



— oo 



for all p, 5 e V'^. In Eq. (281) we introduce 



m 



U{p\ q\ P') = log,, i = 1, 2, 3, (284) 

and combine the integrals from — oo to and from to oo, taking the sign- functions into account. Substituting 
= T we get 

oo 

X f dT r 

aip,q) ^ -CyCAm / — e(Tp^ - 4m^)/i(Tp^ rg^, rP^) 
TT J T I 
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+e(Tg2 - W)/2( V, rg^ rP") + QirP^ - W)/3(V, rg', rP' 



(285) 



Since the anomaly is a polynomial of degree cus + 1 = 2, a(p, g) must be a pure number independent of 
p, q, we can take the limit p"^ ^ and — in Eq. (285), while keeping P^ > 0. Then only the last term 
contributes. Substituting rP^ = s, we obtain 



oo 



a{p,q) = -clcAm / — /3(0,0,s). (286) 

TT J S 



We have for P^ > 4m^: 



which implies 



. . 9 9 9^ 2m 1 - v/1 - 4m2/P2 

/3 p2 = 0, g2 = 0, P^ = , log V ^ =. 

' ^ (27r)5p2 ^l + y^l-4mVP2 

oo ^ 



Am 2 f ds 
a{p,q)=Zj^CyCA J ^ log 



4m2 

Substituting a; = Ath? /s, we get 



1 + V^l -4mVs' 



1 

2i 2 7 1 1 ~ Vl ~ 3; 



1 + 



which shows the mass independence of the anomaly. The further substitution ^Jl — x = z makes the integral 
elementary and we get 

2i 

aij>,(l) = --j^^CyCA. (287) 
Summing up, according to Eq. (282) the axial anomaly for the triangle graphs is equal to 

= __A_c^cA^^A^^^=^"Vg^. (288) 

We have still to investigate whether there exist other splitting solutions which do not have an anomaly 
while preserving all desired properties of the theory. i'^^^^^{p, q) is a pseudotensor of rank two. The lowest 
order normalization polynomial with this property is ~ e^^^^^^^Paqp- But this has already a; = 2, in contrast 
to cOt^ — 0. Hence, renormalization of the VVU triangle does not help. There seems to be a better chance with 
£/iiM2/i3(^p^ g), which is a pseudotensor of rank three with ub — 1. The most general normalization polynomial 
which preserves unitarity and the symmetry in the two V vertices is now given by 

p«^2M3(p^ g) = C£'^i^^'^3a(^^ _ 

where C is a real constant. But this would destroy vector gauge invariance 

p^^P^^^'^^^^^(p,q) = _Ce'^i^2M3«p^^^^ ^ 0, 

which we do not allow for. That means that the axial anomaly cannot be removed by renormalization. We have 
to live with it. In the electroweak theory the anomalies cancel by compensation between leptons and quarks 
([6], sect. 4.9). 

The analysis of the axial anomaly presented above may appear technical at first sight for a reader which is 
not yet familiar with the causal method. However, working with divergent Feynman integrals has some ad hoc 
character, and since the causal method works without divergent, i.e. ill-defined expressions, the computation 
of the axial anomaly presented above is more rigorous than in other approaches and serves as an unambiguous 
consistency check. 
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5.2 Schwinger Model 

The Schwinger model [15, 16] is a popular laboratory for quantum field theoretical methods. As a soluble 
quantum field theoretical model, its nonperturbative properties and relations to confinement [17, 18] have 
always been of greatest interest. It is also possible to discuss the model perturbatively in a straightforward 
way. The interesting features of the model, originally designed to describe QED with massless fermions in 
1+1 -dimensional spacetime, are related to the fact that the massless fermions and the photon field actually 
disappear from the physical spectrum, whereas a "physical" massive scalar field appears with the so-called 
Schwinger mass = e^/ir, which depends on the coupling constant e. For a full discussion of the model 
we refer to the literature (see [19] and references therein). In the following, we focus on the calculation of 
the vacuum polarization (VP) diagram at second order, where the appearance of a mass term in the Schwinger 
model can be traced in the photon propagator by resummation. It turns out that the correct treatment of the VP 
is a delicate task, where the careful discussion of the scaling behavior of distributions given above becomes 
very useful. 

5.2.1 The Causal Approach 

As previously discussed, the S'-matrix is constructed inductively order by order as an operator valued functional 
in the 1+1 -dimensional case 

S{g) = 1 + / d^xi...d^XnTn{xi,...Xn)g{xi)...g{xn), (289) 

n=l 

where g{x) is a tempered test function that switches the interaction. The first order interaction term for QED 
given in terms of asymptotic free fields is 

Ti{x) = ie : ^{x)^^'^{x) : A^{x). (290) 

We note here that the so-called adiabatic limit g{x) 1 has been shown to exist in purely massive theories 
at each order of the perturbative expansion of the S'-matrix [10]. We therefore keep a mass term for the 
fermion fields in the following, and consider the limit m ^ when appropriate. We further note that, of 
course, the properties of 1+1 -dimensional fermion fields do not have much in common with the corresponding 
counterparts in 3+1 -dimensional spacetime, both from a physical and mathematical point of view. 

The interesting second order distribution T2 is constructed by first considering the causal distribution 
D2(x,y) 

D2{x,y)^[T,{x),T,{y)] , (291) 
suppD2^{{x-y)\{x-yf>0} , (292) 
which has causal support. Then D2 is split into a retarded and an advanced part D2 — R2 — A2, with 

suppR2 = {{x-y)\{x-yf>Q,{x'-y')>0} , (293) 

suppA2 = {(x-y)|(x-y)2>0, -(x°-y°)>0} . (294) 

Finally T2 is given by 

T2{x,y) = R2{x,y) + T,{y)T,{x) = A2{x,y) - T^{x)T^{y) . (295) 

For the massive Schwinger model with fermion mass m, the part in the Wick ordered distribution D2 
corresponding to VP 

D2{x, y) = e'idt^^'ix - y) - d^iy - x)] : A,{x)A,{y) : +... (296) 
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then becomes after a short calculation 



d^-^k)= (g.,-^)— , ^ sgn(fc°)9(fc^-4m^) . (297) 

Obviously, df^^ has a naive power counting degree Up — —2 [20]. But the singular order of the distribution 
is a; = [21]. Indeed, applying the definitions from sect. 2.1 to d^'^{k), we obtain the quasi-asymptotics 

limdf = ^ (^g^'^k^ - A;'^r)5(A;2)sgn(A;°) , (298) 

and we have p{6) = 1, hence a; = 0. The quasi-asymptotics differs from the naively expected formal result 



j^^-^, (299) 

which would be ill-defined as a distribution in 2 dimensions. Note that the ^f'^^'-term in Eq. (298) does not 
contribute to the quasi-asymptotics. The reason for the result Eq. (298) can be explained by the existence of a 
sum rule [22] 

DO 

S^m'^ 1 



/ 



SO that the l.h.s. of Eq. (298) is weakly convergent to the r.h.s. In spite of sga(k^), the r.h.s. of Eq. (298) is a 

well-defined tempered distribution due to the factor {g^^k"^ — k^k"). 

This has the following consequence: The (Fourier transformed) retarded part r'^" of d'^'^ would be given in 
the case a; < by the unsubtracted splitting formula 



oo 

{k)^^ [ — - — d^'itk) 



= — U^'^-— - log L= , fc^>4m^A;°>0. (301) 

This distribution will vanish in the limit m — 0, and the photon would remain massless. But since we have 
a; = 0, we must use the subtracted splitting formula 

00 
—00 

= ^ f'^ log — — , k^ > 4:m^,k^ > 0. (302) 

The new local term survives in the limit 777 ^ 0. After a consistent resummation of the second order VP 
diagrams, which we will not discuss here, the well-known Schwinger mass term 777^ = b^/tt appears. Con- 
sequently, the difference between simple power-counting and the correct determination of the singular order 
is by no means a mathematical detail, it is crucial for the proper description of the dynamics of the model. 
The singular order = of the distribution di^'^ implies that a local renormalization is admissible, and even 
necessary to preserve the gauge structure of the theory. 
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An further important property of Eq. (302) is its behavior for k"^ 0. From fg'^ one obtains the corre- 
sponding VP amplitude tg'^ by the replacement m? ^ m? — iO for arbitrary k. It is then straightforward to 
show that the logarithmic term in Eq. (302) behaves for > fc^ ^ like 

log ^ , = log ^ , ~ W-k^/m^, (303) 

%/r^^4^^^-l 1 - ^-kyim^ 

where fc^ < and log(l + z) = z for \z\ <^ 1 was used and the case < /c^ — > behaves accordingly. 
Therefore, we obtain the VP amplitude 

i7{k)=[9^^ -^^)Uk), (304) 

where £2 = ^2,^1'' with t2(^) for /c^ — > 0. Of course, this observation is a direct consequence of the central 
splitting solution used in Eq. (302). 

The causal method provides the most unambiguous guide to the construction of every order of the pertur- 
bative ^-matrix. 

5.2.2 Dimensional Regularization 

We start with the traditional expression for the VP in 1+1 dimensions given by 

where rn? is used synonymously for rn? — iO. According to the recipes of the dimensional regularization 
procedure, we consider the trace of Eq. (305) in n dimensions 

Z^^(l,^ _ _on/2 f d> (2-n)(p^-p/c)+nm^ 

J^W- y (27r)-(p2-m2)[(p-A;)2-m]- ^ ^ 

We will now proceed in two different ways. First, we perform a naive dimensional regularization by taking the 
limit (2 — n) — > in Eq. (306) before performing the integral. Then only the term 

_ 2-/2, / ^> GOV) 

remains. Inserting the Feynman parameter integral 



1 

1 f dx 



ah J [ax + 6(1 — x)\^ 


and using the 't Hooft-Veltman formula [13] 

^ r dy ^ , , /, r(a-n/2) 

°' J (p2 - 2pA; - m2)« T{a) 



(308) 



T{a) (A;2 + m2)«-V2 
for a = 2 and n = 2, we obtain 



(309) 



^°^=- .(i (310) 



and hence 



1 



I 



m^dx 



^m(^) = - / -T^ ur2 2- (311) 

TT J x[l — xjk'^ — 
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The important observation is that obviously 



^ (^) ^ -V^r for ^ ^ 0, (312) 

i.e. we are left with the same problem as in Eq. (301) that Eq. (312) reproduces only "half the VP amplitude. 
Note that the result differs by a factor 27r from the result derived in our causal framework, since there a 
symmetric definition of the (inverse) Fourier transform has been used. 

We now renormalize Eq. (305) properly according to the prescriptions of 't Hooft and Veltmann and 
include all terms in the integral. The full tensor structure is given by 

1 

V-l'^; ^ J {2t:Y J [p2_2pA;(l-x) + A;2(l-x)-m2]2 " ^^^^^ 



To perform the integrals, we use the 't Hooft- Veltman integrals 



I 



(p2 - 2pk - rri?Y 



hk^", (314) 



for a = 2 and n — > 2 and obtain 



1 

V(^) = -^^^ j x{i- x)k^ -m^ ^^^^^'' - (1 - x)k'^k^{g^aguf3 - g„vgafi + gnngva) 



+[(1 - xfk'^k^ + 7^{x{l - x)e - m^)g''^]{g^^g,fs - g^.gap + g,.pgua)]- (316) 

Interestingly, the integral above is finite due to a cancellation of the dimensional pole (2 — n)~^ by a factor 
generated by 

gmQi^P - QnvgaH + gupgva = (2 - n)g^^, (317) 
since gj^ = n. This way we arrive at the gauge invariant amplitude 

i^,u{k) = {k'^g^, - k^K)t{k), (318) 
x{l-x)dx ^3^^^ 



TTJ X[l 



x)k'^ — vn? 



with vanishing trace in the limit /c^ ^ 

t>W = rt(.) = -i^ + „(^)^ 0. (320) 

Both results in the dimensional and causal regularization scheme are consistent, however, we observe that the 
regularization of distributions must be performed with due care. 



38 



5.3 Scalar QED in 2+1 Dimensions 



In this section we illustrate how gauge invariance is automatically preserved by dimensional regularization 
by using scalar quantum electrodynamics (sQED) in one time and two space dimensions as a an example. 
As expected, both the causal method and dimensional regularization lead to compatible results, although the 
underlying premises on which the two methods are based and the resulting perturbative description of the 
model theory display a rather different behavior. 

The traditional starting point of any quantum field theory is a Lagrangean containing coupled classical 
fields describing the interaction. After quantization, 5'-matrix elements or Greens functions are constructed 
with the help of, e.g., Feynman rules. One should accept that this point of view is obsolete to some extent 
within the causal approach, which is based on a description of the ^-matrix by the help of free fields and well- 
defined interaction terms which can be expressed as a sum of normal-ordered products of free fields. However, 
the Lagrangean provides a formal tool to express the classical structure of the theory, but one should keep in 
view that physical theories are always subject to quantization. 

The scalar sQED2_,_i Lagrangean is given by 

C = -^F^uF'"' + (a^ - ieAi,)^\d'' + ieA^)ip - mVV, (321) 

where describes the electromagnetic field with the field strength tensor Fni, — d^Ai, — d^A^j, and </? a 
charged scalar meson field with mass m and electric charge e. The coupling constant e has the dimension of 
an energy in three-dimensional spacetime, and consequently SQED2+1 is superrenormalizable by naive power 
counting. 

The Lagrangean can be decomposed according to 

^ ~ ^em + ^matter + ^inti (322) 

where the interaction part Cint is given by the minimal coupling of the electromagnetic current to the electro- 
magnetic potential 

Ant = -Jm^^' = ^^v'^ dn<p- e^(p^ipAu,. (323) 
It is now straightforward to construct the Hamiltonian interaction density 

Hint = ie<p^ <pA'' - e'^ip^ipA^Ai' + eVV(>l°)'- (324) 

Obviously, this expression is not manifestly covariant. It has been shown in [23] that the non-covariant term 
—e'^^p'^^p{A^Y is canceled in the full perturbative quantum field theory by a local normaUzation term appearing 
in the so-called seagull graph at second order in the coupling constant which is sesquilinear in the meson and 
bilinear in the photon field. This observation applies to scalar sQED in n + 1 dimensions in general, and a 
detailed discussion of sQED3_|_]^ with respect to gauge invariance can be found in [24, 7]. 

It must also be mentioned that scalar sQED has a pathological infrared behavior both in 3 + 1 and in 2 + 1 
dimensions. E.g., the contraction of the two photon lines in the seagull graph with the two photon lines in 
another seagull graph gives rise to a (/^VV'/^'iiiteraction which is highly singular at short distances, and it is 
generally accepted that this leads to a transmutation of the original underlying perturbative theory. However, 
in this work we will focus on strictly perturbative aspects of the theory, which are well-defined at every order 
of the coupling constant e in our case. 

5.3.1 Causal Approach 

The crucial difference between the two approaches discussed in this work is the following. In the causal 
approach, the interaction Hamiltonian density is given by the normally ordered product of free quantized 
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fields, whereas in the dimensional regularization ansatz, all kinds of UV divergences including terms which 
stem from contractions of fields at the same space-time point are taken into account. 
In the Feynman gauge, the free photon field A^{x) fulfills the wave equation 

DAi,{x) = (325) 

and the translation invariant distributional commutation relations 

[A^{x),A,{y)] = [A^{x - y),A,{0)] = ig^.o'^^Xx - y), (326) 

where 

^S^'H^) = /" d'pS(p')e(p')e-'^' (327) 



applies for photonic contractions without time-ordering. For the scalar field, we have the contractions 

^^\y)^-iD^^\x-y), (328) 
i^^iy) = +tD!;^>ix -y) = +tD^^\y - x) = -tD^^\x - y), (329) 

where 

D^+\x) = — ^ / £p5{p^ - m2)e(p°)e-^f^ (330) 
(27r)^ J 

One may choose as a starting point a Hamiltonian density which is given by the normally ordered products 
of free fields ^ 

Uintix) = -{ie^\x) ^{x) + e^^\x)^{x)A^{x))A^{x) (331) 

and X is an element of 2 -|- 1-dimensional Minkowski space. The perturbative ^'-Matrix is then constructed 
according to the expansion 

S=l + Y. / dxl..dxlT{Hint{x,)Hint{x2) ' .■■ ' Hint{Xn)}, (332) 

n=l ' 

where T is the time-ordering operator. 

In the causal approach, we use only the first order term (in the coupling constant) of Eq. (331) (which is 
motivated by first order interaction term appearing in the corresponding Lagrangian) given by 

Ti(x) = e : ^\x) v{x) : A'^(x). (333) 

Thus, the primed distributions can be written as (taking into account that Ti = — Ti) 

A!^{x,y) = -e^ : (^t(^) (^(^) ^(y) A''{x)A''{y) (334) 

R'^ix, y) = -e^ : (^t(^) ^^y) ,, (^t(^) ^^^) ,, Ai'{x)A^{y) :, (335) 

the causal distribution is again 

D2{x,y)^R'^{x,y)-A'^{x,y) = [Ti{x),Ti{y)]. (336) 
As an example, from Wick's theorem we have 

: ^t(^)9V(^) :: ^^{y)d''ip{y) : 

=: ip\x)d^^ix)^^iy)d''^iy) : +^^ix)d''^iy) : id^Di^\x - y) 
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- : d^^{x)^\y) : idlD'^-\y - x) + d^D^+\x - y)d;D^-\y - x), 



with Dm^ the positive/negative frequency part of the Pauli- Jordan distribution. 

We focus now on the vacuum polarization diagram, corresponding at second order to 

D2{x,y) = -e^:A^{x) A,{y) : 
• - x) d^D(-\x -y)- d^D^'^y - x)d^,D^^\x - y) 

+9,^L>W(y - x) d^^D^-\x - y) - d^,D^-\y - x)d;D^^\x - y) 
-d^d-^D^r^Hy - ^) Di;;\x -y) + d^^dlD^-\y - x) D^^\x - y) 
-D^^\y - x)d^^dlD^-\x -y)-r D^-\y - x)d^^dlD^^\x - y)] . 

For further calculations, we get rid of external fields and change to momentum space 



D2{x, y) =: - y) : A,{x)A,{y) : with d^^{p) =: -r U 



p2 



d{p). 



(337) 



(338) 
(339) 
(340) 
(341) 
(342) 



(343) 



In a straightforward manner, one derives (p = 



P 



d{p) = i^e(p - 4m^)sgn(po) 1 



Arm? 
p2 



(344) 



Obviously, this distribution has singular (and power counting) degree u — 1. 

The retarded distribution f(p) follows directly from the splitting formula (first for p = 0), — Pq > 0) 



iO)^+i (po -ko + iQ) 2 



r(po) = ^(27^) p%^ J dko 
This dispersion integral can be written 

oc 

f(po) = 2i7r^Po j 



1 - 



(345) 



A;2(p0 - fco + iO) A;2(p0 + fco + iO) 

oo 



4m" 



2i7r^Po 



+ ^0)2 - 



oo 

/ 

4m2 



ds 



Am? 



g3/2 p2 _ g _|_ ^QpO ' 



(346) 



where a substitution s = kg was used. For the C-number part r'2 of i?2 follows 

r'iPo) = -|e(p^ - 4m^)e(-po) (l - ^) 



and therefore 



00 

iipo) = ^(Po) - r'ipo) = 2i7rVo ^ 

4m2 



(is s — 4m2 

g3/2 p2 _ _g _|_ 



(347) 



(348) 



The integral can be evaluated by standard methods, and proper analytic continuation of f(po) leads to the result 



p2 



2m p2 — 4m2 2m + p 



p 



+ 



2p2 



log 



2m — p 



(349) 
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5.3.2 Dimensional Regularization 

Applying the well-known Feynman rules to the vacuum polarization diagram leads to 



3 ■2h-i'p" + ■2h"pi' - jfjf - Ikt'k" 



f"" = / d'k ' \ , , (350) 

J (A;2 - m2 + iO) [{k - pf - m2 + iO] 

using the definition of the Feynman propagator 

(0| T y{x)^\y)\ |0) = -iD^{x - y) . (351) 
The integral can be split into a scalar, vector and tensor part: 

t^'' = 2p'^X^(p) + 2p%^(p) - p<^p'I,{p) - 4Xr(p) , (352) 



with the definitions 



1 




-TV? + iO) [{k - 
k^" 


pf 


— m? + iO] 




-w? + iO) [(A; - 
ki'k'' 


pf 


— m? + iO] 




-w? + iO) [{k - 




— + iO] 



(353) 
(354) 

(P - m2 + iO) [(A; - p)' - m2 + iO] ■ 

We first calculate Xi. Using the Feynman parametrization 

1 

= j da [aA + {1- a)B] 



and performing the momentum translation k^ k^^ + ap^, the integral can be written as 

1 

I,(p)^/A/.„[*^-,„^ + „(l-ay+«]-. (357) 



Changing the integration dimension to D = 3 — 2e and using the general relation 



/ 



\m-r- 2 



(356) 



d k——^ = TT 2 — ^ — ^ ^ , (358) 



the momentum integral can be carried out (with the trivial limit e — > 0), after which the cc-integral becomes 
simply 

1 

Iiip) = ITT / == = — log , (359) 

J y/a{a-l)p^ + P 2m - p 

where p = a/p^. 

Note that Eq. (359) is indeed valid for arbitrary p, when is substituted by — iO. For > Am'^, 
the logarithmic term contains both a real and an imaginary part, for < < Am'^, the logarithm is real. For 
space-Uke momenta p^ < 0, the logarithm becomes purely imaginary, but also the prefactor 1/p — 1/ — 
—i I —p^, since the integrand in Eq. (359) is real in this case. For p^ < 0, one can also write 

2:1 (p) = log = ^=arcsm . „ , „ , (360) 



'p 



2m - ^ I 1 - 4mVp2 / ' 
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but we will maintain the shorthand used in Eq. (359) in the following. 
The same procedure as above can be applied to I2, leading to 



1 



(361) 



Integrating symmetrically makes the integral proportional to disappear, leaving 



1 



da 



m'^p^ 2m + p 

log; 



2p 



2m — p 



(362) 



-^3 



T3 Finally, we use Feynman parametrization and momentum translation invariance in order to obtain for 

1 

X^^(p) = J d^k J da [k'^k'' + a {k'^p'' + p'^k'') + aVp^] 



X [k^ -m^ + a{l - a)p^ + iO] ^ , 



(363) 



where the integrals proportional to an odd power of k vanish. 

Rewriting k'^k^ — g^'^k'^/D and performing D — 3 — 2e dimensional integration, we arrive at a finite integral 
for e ^ 



ir{p)^in'l 



da 



g^^'^aia - \)p^ + + p^p"" 



a 



a[a — + m^ 



which can be evaluated in a straightforward manner to give 



-p 



p2 



2m p^ — 4m^ 2m + p 



P 



+ 



2p2 



log 



2m — p 



(364) 



,„Am p^p" , 2m + p 
p 2m — p 

Now, all integrals can be combined in the manner expressed by Eq. (352) to give the full VP 



p2 



2m p^ — 4m^ 2m + p 



p 



+ 



2p2 



log 



2m — p 



+Ai'K'^g'^''m . 



(365) 



The interesting observation is given by the fact that the result obtained so far does not exactly match the 
gauge invariant result derived in the framework of causal perturbation theory. In fact, we have to include the 
one-loop contracted seagull graph displayed in Fig. (4), which contributes to the photon-photon transition 
amplitude like the VP diagram displayed in Fig. (3) as well. 

Calculating the formal contribution of the fermion-line self-contraction of the seagull graph gives (with the 
correct normalization factor) by the help of Eq. (358) 



seagull 



^-2g 



d^k 



^2 _ ^2 _|_ 



(366) 



which exactly cancels the local term in the scattering matrix element Eq. (365) of the vacuum polarization, 
which only appears if one uses dimensional regularization. 
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Figure 3: Vacuum polarization 



Figure 4: One-loop contracted seagull graph 



5.3.3 Pauli-Villars Regularization 

Xi, according to naive power counting, is convergent and can be calculated using Feynman parameters 

l^{p) = — log ^. (367) 

p 2m — p 

For I^, we change the propagators according to the procedure 

' ' ' (368) 



q2 _ ^2 _j_ q2 _ ^2 _|_ g2 _ ^2 _|_ ' 

leading to 

2 J (p_^2)(p_^2)p_p)2_^2][(^_p)2_/^2] ■ 

Subtraction at p = gives 

IJ(p)=IJ(0)+Jnp), 
X,^(0)=/a. -A')' ^ 



{k^ - m2) (/t2 - A2) 

+ r.,, „o!r.,, .0 .0. !^ • (370) 



- m2) (P - A2) [(A;-p)2_m2]p_p)2_A2] 

Taking the limit A ^ 00 yields 

I^{0) = 



I^ip) = J d'k 
Using Feynman parametrization 



{k"^ — m?) [(A; — p)2 — m2] 



1 a 



A^J = J da J P)A + PB + {1- a)C]-^ (371) 
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and momentum translation leads to 



1 a 

= J £k 2 J da J dp 







x[2p-kk^' -{1- 2(5)p^k'' + 2/?p> • A; - (1 - 2/3)/3pV] 

X [k^ - - 1)/ - + iO] . 
Again, the integrals proportional to an odd power of k disappear and we find 

"2 



1 a 



J da J dp J 



d-'k 



-k' - (1 - 2p)pp' 







which can now be Wick rotated and integrated over k to 

2 



1 a 



(372) 



(373) 



= m'^p^ I da 



a 



•y/ a{a — l)p^ + w? 



i-K^p^ , 2m + p 

— log = 

2p 2m — p 

The same procedures may be applied to Xg"^, too. One has 



X^ip)- 



-^3 



-P 



r- 



p2 



2m p^ — 4m^ , 2m + p 

— + — 7r-2 — o 

p 2p^ 2m — p 



+ 



p^p^ 2m + p 



2m — p 



Combining all three integrals obtained by Pauli-Villars regularization, one obtains 



p. 



2m p^ — 4m^ , 2m + p 

— + — — o 

p 2p^ 2m — p 



(374) 



(375) 



(376) 



which differs from the solution with dimensional regularization by a term Ain'^g'^'^m, i.e. gauge invariance is 
preserved "by hand" by proper normalization of the photon-photon transition amplitude in the present case. 
However, this is only true if one starts with a normally ordered first order coupling and neglects the one-loop 
contracted seagull graph from the beginning. 

Again, we observe that the calculations differ strongly in different approaches. In the causal method, the 
non-trivial part of one-loop calculations is a one-dimensional, finite dispersion integral. As a general remark, 
we point out that the causal method, as well as dimensional regularization, has some particular advantages 
for gauge theories because it does not use a cutoff which breaks gauge invariance, an observation which is 
illustrated on a rather basic level in this sect. 5.3.3. In a more general setting, it works in a fixed number 
of physical dimensions so that problems originating from axial couplings, which are also related to the '75- 
problem' in dimensional regularization. 
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Figure 5: Lowest order scattering Figure 6: Bremsstrahlung diagram 

6 Infrared Divergences 

The infrared structure in the causal approach differs strongly from other approaches. Whereas infrared di- 
vergences show up as poles in dimensional regularization or as divergences in the artificial mass parameter 
introduced for the originally massless fields contained in the theory under consideration (a strategy commonly 
used e.g. in the Pauli-Villars approach), they are automatically regularized by the test function g in Eq. (189). 
From a mathematical point of view, this is the most natural formulation of the infrared problem, since the T„'s 
are operator- valued distributions, and therefore must be smeared out by test functions in <S(E^"), the Schwartz 
space of functions of rapid decrease. The test function g E 5(R^) plays the role of an "adiabatic switching" 
and provides a cutoff in the long-range part of the interaction, which can be considered as a natural infrared 
regulator. An appropriate adiabatic limit g ^ 1 must be performed at the end of actual calculations in the right 
quantities (like cross sections) where this limit exists. Strictly speaking, the perturbative S-matrix according 
to Eq. (5) does not exist for (7 = 1 for many theories involving massless fields. This observation is also closely 
related to the notion of infraparticles [26, 27]. The existence of the adiabatic limit for physical observables 
is a non-trivial issue of a theory and is not automatically guaranteed at every order of the theory from the 
mere lowest order definition of the interaction. Below, we will prove the existence of the adiabatic limit for a 
scattering cross section at fourth order in g for a model theory. 

Note that introducing a mass as infrared regularizator for massless fields is a questionable procedure, since 
it is unclear whether the original massless theory is restored by taking the massless limit of the massive theory, 
which may suffer from potential problems like, e.g., broken gauge invariance. 

In order to demonstrate the causal approach to the infrared problem we consider a theory in 3-1-1 spacetime 
dimensions, called totally scalar QED in the following, where a massive scalar charged field is coupled to a 
massless scalar field, in close analogy to Eq. (216). The corresponding scalar particles will be called meson 
and photon in the following. The theory is defined by the first order coupling term 

Ti(x) = -ie : ip\x)'^{x) : \Aq{x) + A"'^*(x)] = -ie : '^\x)^{x) : A{x), (ill) 

where A^^^{x) denotes an external C-number field and Aq{x) the quantized massless neutral scalar photon 
field. For dimensional reasons, the coupling constant e has the dimension of an energy or an inverse length. 

In the following, we consider the scattering process of the meson off the external field, according to Fig. 
(5). At first order in the coupling constant, the matrix element for the scattering of a meson with mass m and 
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initial state momentum q and a different final state momentum p is given by 



M;) = (/|5«|i) = (0|a(p-)5«at(5-)|0), 



(378) 



where S^^'^ = ~ie J d'^xi : ip"^ {xi)ip{xi) : A^^*{xi) and \i) = a){q)\Q) and |/) = a){j))\{)), i.e. a and a) denote 
the annihilation and creation operators for one charge type of particles of the meson field, and the vacuum shall 
be denoted by |0) in this section. One obtains {k^ — {k"^ + m^)^/^ etc) 



ax^ / ^«a«_g-ifexigifc'xi^Q|^(^ . Q(^)«t(^/) . a\q)\Q)A''^\x^). (379) 
2y /co^o 



(27r) 



Exploiting the commutation relations and the distributional identity J dP'xe^^^ — (27r)"5(A;) leads to 



/I r rl''>hrl''>h' ^ ^ 



—le 



(27r)3y ^^^2vp^ig;^^ ^ ^ ^^^^ 

Introducing the Fourier transform of A'^^^{xi), we have 



—le 



2(27r)3^^ 



dSie-*''^ie*^'^M^"*(xi) (380) 



—le 



— ie 1 

2(27r)3v»(2^ 

— ie 1 
2(27r)3^/M^(2^ 

and the first order matrix element becomes 



c^^xie-^^^ie^^'"^^^ y d^k"e-'^"''^A'''\k") 



d^xi j d^k"e-''^''^e'^^e-'^"''^A^''\k") 
d^k"{2T:)%-q +p- k")A'^\k"), 



—le 



A^^\p - q) 



2(27r)v» 

We assume for the moment that A^^^{x) is a Coulomb potential 



(381) 



(382) 



\x\ 



(383) 



For the sake of completeness, we calculate the corresponding cross section in detail. In the present case, we 
have 



(1) ^ le 



1 



{2n) 



d x- 



J{p-q)x' 



\X\ 



le 



(27r)32VP^ 



-i{p-q)x 



\X\ 



(384) 



The space integral can be evaluated as follows 



/ d^'xe-'^^-^^^^ = -jJ-^ [ d^x^Ae 

J Fi [p-Q) J m 



-i{p-q)x 



d^x 



(4)' 



-i(p-q)x 



(P- q)' 



d^x(^ - 47rS^^\x)y-'^P-^' 



An 



(P- ^)^' 



and thus 



c(i) ^ _ 



{27tS{po - qo)) 



{p-q}^' 



(385) 
(386) 
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The transition rate dR from the initial state to a final state within an infinitesimal phase space volume d^p is 
given by (T denotes a large time interval) 

Of course, the square of the ^-distribution above is ill-defined, since one should work with wave packets in 
order to get well-defined expressions. However, for the moment we content ourselves with Fermi's trick and 
perform sovae, formal manipulations, starting from 

T 
2 

27r5(po - ?o) = lim / dte'^°-'^''^K (388) 

T^oo J 



The intuitive argument is that this expression is non-zero for po = so that one may replace one 5-distribution 
by 

T 
2 

27r5(po - go) = lim / dt (389) 

T^oo J 



T 
' 1 



and for large T one has 27r5(po — ^o) = T- Then dR becomes 

2 _/3 

dR = — ^ {2nS{po - qo))j^\,. (390) 

(27r)%go (P - g)^ 

The cross section is given by the ratio of the transition rate dR and the flux of incoming particles, given by the 
initial state expectation value of the operator 

f^l: ip\x)df'ip{x) - ip{x)d^'ip\x) :, (391) 



leading to 

n = 

for the chosen initial state. The differential cross section is therefore 

d(j = ^ = ^—4 (27r5(po - q^)MdW^ 



(392) 



e 



2 



Since poc?Po = \p\d\p\ and |p| = we obtain 

o2 



<^(po-go)|pTc?|p1(^^^. (393) 



da = j^^-^S{po - qo)dpodn (394) 



and finally, performing the integral over po, 

da f ^, , , 



i.e. we basically recover the Rutherford cross section. 



(395) 



P0=90 
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6.1 Bremsstrahlung 

We now consider the case where the scattered meson emits a soft photon with four-momentum k according to 
Fig. (6). The scattering matrix at second order is given by 



[-lef 



d xid X2 : (p{xi)(p\xi)A{xi) :: (p{x2)(p\X2)A{x2) : g{xi)g{x2) 



(396) 



where the contraction symbol denotes one possible Wick contraction of massive fields. Since one has two 
possibilities to contract the massive fields, the relevant bremsstrahlung term becomes 

^(2) = (_ie)2 j d^x^d^X2{-i)D^{x^ - X2) : v{x^)A{x{) :: ip\x2)A{x2) : g{x{)g{x2) + . . . (397) 

and the dots denote other terms from the Wick ordering of jS"^^^ which are irrelevant for the present case. The 
external field operators have the form 



^\X2) 



1 P d^ k" 



A{X2 



and we obtain the bremsstrahlung matrix element S^^'^^ — (0|a(p)c(/c)5'a^(g)|0) after some calculation 

1 



(398) 

(399) 
(400) 
(401) 



gWems ^ (_ •g)2 j d''x^d^X2{-i)D'^{x^ - X2 



+ 



^31,(4) j3 7,./ ^3^,/// ^ ^ ^ ^ ■ 



(402) 



(27r)9/2 



d xid X2D^{xi - X2)g{xi)g{x2) 



VSqoPoko 



A^"*(X2) + 



y/8qoPoko 



Inserting the Fourier transforms of g'(xi), g'(x2), A^^*(x2) 



(403) 



d%e-''''-'g{k,), g{x2) 



(27r) 



d^kie-'^''''A^^\ki), A^''\x2) = 



(27r) 



^4^^g-ifc2X2^ea;t(^^) (404) 
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and the Fourier transform of the Feynman propagator leads to 



nbrems 



-le 



(27r)W2^8^^(27r)4(27r)^ 



d^xid'^X2 



Q—ik"x\ Q—iqx\ ^ikxi ^—ikixi ^ik"x2 ^ipx2 ^—ik'x2^—ik2X2 

„—ik"xi „—iqx\ „—ik'xi „—ikixi Ak" X2 c,ipx2 „ikx2o—ik2X2 
e e e e e e e e ■extnJ\ 



^ext^j^ 



g{ki)g{k2). 



(405) 



Note that we sometimes omit the iO-term of the Feynman propagators for the sake of brevity. We evaluate all 
trivial integrals and arrive at 



nbrems 



2(27r)i3/2V2goPo^o 



d'^kid'^k2 



1 



1 



{q — k -\- kiY — m? {p + k — k2)^ — m? 



A'^\p-q + k-kr-k2)g{k^)g{k2). 



(406) 



ki can be replaced by k2 



obrems 



-le 



2(27r)i3/2^2^^ 



d'^kid^k2 



1 



1 



{q — k + k2y — m? {p + k — ^2)^ — 



A'-\p-q + k-k,-k2)g{k,)g{k2). 
Now we investigate the adiabatic limit by first replacing g{ki)g{k2)hy 

1 . fki\ 1 . fk2\ 



(407) 



(408) 



corresponding to the replacement g{x) — > go{ex) in real space, and taking the limit e — > 0. go{x) is a fixed test 
function in <S(R^) with 510(0) = 1, so that go{ex) ^ 1 for e ^ 0. Note, however, that 1 ^ »S(R^). Thus we 
have 



obrems ^ ^^^^ f ^ ^1 



d^k2 



2(27r)13/2^/g^^ 
1 



(p + A; — A;2)^ — 



(g — + ^2)^ — 

^"1 \ . / k2 



or 



obrems 



+ 



2(27r)i3/2^g^^ 

\ext 



A^'^\p-q + k-k,-k2)go{j)go{j) 

= I d%d^k2 
J 



(409) 



(g — A; + eA;2)^ — 



A^'^^p - g + A; - e/ci - ek2)go(ki)go(k2). 



(410) 



(p + /c — eA;2)^ — 

Envisaging the limit e ^ 0, we can neglect the e-dependent term in the argument of A'^^^ and perform one 
trivial integral. We obtain the result 

1 1 



obrems c(l) 

~ (27r)V2^ 



2p{k-ek2) 2q{k-ek2) 



goih 



(411) 



where we used the first order matrix element S^^^. The cross section follows, using the lowest order cross 



section 



da{k) da^^^ 



dn dn (27r)72|A;| 



d'^h 



1 



1 



-2q{k-ek2)) 2p{k-ek2), 
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^0(^2) 



dn 4(27r)72|A;| 



1 



1 



1 



{q{k - eki)){q{k - e/cg)) 
1 



+ 



For well-known physical reasons, one has to integrate this cross section over the photon momenta up to a cutoff 
ojq, assuming that photons with momenta < ojq are not measured. The four integrands can be rewritten by the 
help of the Feynman trick, leading to 



da 

dn 



d-^k 



|fe|<a;o 



4{27iy2\k\ 



d^kid'^k2 



dx 



1 



1 

/ 



dx- 



1 



[q{k — eki)x + q{k — ek2){\- — x)Y J \p{k — eki)x + q{k — ek2){\- — x)]^ 





— I dx 



[q{k — eki)x +p{k — ek2){l — x)Y 



+ / dx 



[p{k — eki)x +p{k — ek2){l — x)]'' 



9o{ki)go{k2). 

(413) 



We first consider the fist part of the integral above, the last term can be treated in an analogous manner. 



Ii — dx d kid k2 



d^k 



|fe|<a;o 



= dx d kid k 



2\k\ [q{k - eki)x + q{k - ek2){l - x)f 

\k\d\k\{-d cos e)d(f) 



|fc|<<^o 



2[go|A;| - cos^ - (qeki + (1 - x)qek2)] 



9o{ki)go{k2) 



-9o{ki)9o{k2) 



— I dx I d'^kid'^k2 



\k\d\k\ 



-1 



27rdcos^ 

[-|^|^| cos^ + qo\k\ - {qeki + (1 - x)q€k2)] 



-9o{ki)9o{k2) 



1 UJO 

^2tt J dx J d%d% J 



d\k\ 

m 



1 



(Qo + \Q\)\k\ - (qeki + (1 - x)qek2) 



ilo - \Q\)\k\ - {qeki + (1 - x)qek2) 
9o{ki)go{k2) 



27r I dx I d^kid% 



2k1 



log 



— * 

(go - \Q\)\k\ - {qeki + (1 - x)qek2) 



0)0 



log 



{Qo - \^)\k\ - {qeki + (1 - x)qek2) 



9o{ki)go{k2). 



(414) 



{% + |?1 

Focusing on the divergent part in the expression above, i.e. examining the logarithms in Eq. (414) for \k\ — 
only results in 



I[^2tt J dx J d^kid^k2^^ 



(go + |g1) (?o-k1) 
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log 



{qeki + (1 - x)qek2) go{ki)go{k2) 



2% dx d kid k2 



log|e| + 0(l) go{ki)go{k2) 



(2 



TX 



= -27r dx d kid k2^log\e\go{ki)go{k2) = ^log|e|, 



(415) 



such that we obtain the divergent contributions to the cross section from the first and the fourth part in the 
integral Eq. (413) 



Ii + l4 = S^log e 



(416) 



Finally we calculate the second integral appearing in Eq. (413), the third integral can be calculated analogously. 
We have 



h^- dx d'^kid'^k 



d^k 



|fe|<a;o 



2\k\ [p{k - eki)x + q{k - ek2){l - x)] 



r^o (^1)50(^2) 



1 

^- J dx J d%d% j 



d%d^k2 

° \k\<uio 

1 OJQ 1 

•I;,, j4i 



-dcos 



= —271 dx d kid k2 



\k\d\k\ 



2[p{k — eki)x + q{k — eA;2)(l — x)] 



dcos9 



^go{ki)go{k2) 











2 J [{x{p — q) + q)k — {xpeki + {1 — x)qek2)]'^ 
-1 

We abbreviate x{p — q) + qhy Q below. 



goiki)goik2). (417) 



1 1^0 
'4l j4 



I2 = -2tt / dx / dUiciUs 



\k\d\k\ 



dcos9 







2 i [go|A;| - |g||A;|cos^-(xpeA;i + (l-a;)geA;2)]2 



go{ki)go{k2) 



1 



= -2j- / dx I d^k,d% 



1^0 
" d\k\ 



1 



[{Qo + \Q\)\k\ - {xpeki + (1 - x)qek2)] 



[{Qo - \Qm - {xpeki + (1 - x)qek2)] 
go{ki)go{k2) 



= -27r 



d%d% 



log 



2|g| 



(go-igi) 



log 



— * — * 

{Qo - \Q\)\k\ - {xpeki + (1 - x)qek2) 



{Qo - \Q\)\k\ - {xpeki + (1 - x)qek2) 



{Qo + m 

Again, we evaluate only the infrared divergent part of the expression above. 

1 



go{ki)go{k2). (418) 



I' = -2n / dx / d^kid^k 



2|Q| 



(Qo + igi) (go-igi) 
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log 



{xpeki + {l-x)qek2) go{ki)go{k2). (419) 



Performing the integral over ki und k2 results in 

1 



7^ = 27r 



log|e|+0(l)^o(A:i)^o(A;2) 



= (27r)^log|e| j dx-^ = (27r)Mog|e| J dx- 



x'^ip^ — 2pq + + 2x{p — q)q + q"^] 



Now we use the fact that pq — + ^P^ for P — p — q 



PQ = PoQo -pq^E-pq^m+ ' ' — 



^ ^ p2 

pq — m + 



2 2 



2 ' 



/2 = (27r)Mog|e| y dx 



1 



[_p2^2 + p2^ + ^2]2 



1 r 1 

■(27r)^log lei— / dx ^ — 



We substitute y — x — \ und a — ^ + \ and obtain in a straightforward manner 

1 

2 



(27r)^log|e| 



2m V 4m2 

1^1 

For the sake of convenience, we set 6 := so that 



log 



\P\ 
2m 



p2 



Aw? 



4 = (27r)Mog|e| 



Combining I 2 with the third integral in Eq. (413), we have 

7^ + 7^ = (27r)^log|e| ^ 



log 



=6\/rTF 



m- 



log 



6 + Vl + 62 



and the full infrared divergent part of the bremsstrahlung cross section from I[, I2, und I '4^ is 



dfjbrems ^^(1) g2 



1 + 



log 



log|e|. 



6.2 Self Energy 



(420) 



(421) 
(422) 



(423) 



(424) 



(425) 



(426) 



So far we considered the infrared divergences in the bremsstrahlung cross section, which represents a fourth 
order contribution in the coupling constant e to the inclusive cross section. To see how the logarithmic infrared 
divergences ~ log(e) compensate, we have to investigate the third order Feynman diagrams according to Figs. 
(7) and (8), since they combine with the first order diagram in Fig. (5) according to 



da 

dn 4th order 



~ I'-'filithorder 



\sf^n2Re\S^S^\ 



(427) 
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j^ext 



Figure 7: Self energy diagram 



Figure 8: Vertex diagram 



We consider first the self energy diagram. The relevant contributions to the S-matrix are given by six equivalent 
variants of the contraction 



6 



d xid X2d X3g{xi)g{x2)g{xs) 



I 1 I 1 1 

: ip{xi)Lp\xi)A{xi) :: (p{x2)(p\x2)A(x2) :: ip{x3)Lp\x3)A{x3) : . 

The corresponding third order S-matrix contribution is therefore 

^(3) = (_,e)3 / d^Xid^X2d%{-t)D'^{xi - X2){-t)D'^{x2 - X3){-t)DUx2 - Xs) 



(428) 



: ip{xi)ip\x3) : A{xi)g{xi)gix2)gix3) + • • • , (429) 

the massless Feynman propagator Dp coming from the photon propagator. Calculating formally the matrix 
element with the corresponding initial and final meson states scattering off the external photon field leads to 

"3 ^ A''''\p — q — ki — k2 — ks) 



Qself 



2(27r)37g^ 



d'^kid'^k2d'^k3 



1 



(27r)2((j9-fc2-A;3)2-m2) 



g{ki)g{k2)g{k3) 



(27r)6 k^{{k - {ks - p)y - ^2) J ^v-^^v-^^v-^- ("^30) 
However, the integral above contains a UV divergent part, which must be handled properly. The formal integral 

1 



S'^-(p2,A;3,e)= / d^k 



(431) 



must be regularized or treated within the causal framework. The finite result for the self energy diagram is [28] 



1 H ^ log 



C 



(432) 
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where C is a free normalization constant, and we must replace the formal integral above by 
In the following, we replace again g by and perform all trivial integrals 



Sbrems^g(l)_J^ / d%d% ^ 



ITT 



The integral over k2 gives 



(27r)2 

.2 



1 



(27r)2 -2pe(A;2 + A;3) 



9o{h)go{h). (434) 



{2tiY J ""-2peh 
Now we replace C by (f^ (C - 1) 

4— = / d%. 



ITT 



1 + 



2pek3 r f2pek. 



log['^)-ine{-2peks)]]+C' 



go{ks). (435) 



1 iTT^ 



S 



2pek^{2TiYi 

1 71^ \2pek 



log - z7re(-2peA;3)l + c]Uh) 



(1) / ^4. 

(27r)2y '2pefc3(27r) 



log|e|+0(l) +C goiks) 



MP ./l.. „ log|6| / C^%^o(^3) + 4. 



e2 



(436) 



The bremsstrahlung diagram contained only logarithmic infrared divergences. For this reason, we choose 
C = in order to avoid a 1/ e-divergence in the self energy diagram and obtain 



2 1 

qbrems _ cW 1,1 

-^f^ (27r)2 4m2^°S'''' 
and the corresponding contribution to the cross section follows from 

% ^2Re\Sf^Sf^\ 

dQ 4th order ^ ^ 



(437) 



and is given by 



6.3 Vertex Function 



flfjbrems ^^(1) g2 



dQ,div dO, 2(27r)^m' 



log|e|. 



(438) 



(439) 



Finally, we consider the vertex diagram according to Fig. (8). Formally, one obtains the expression containing 
the UV divergent scalar vertex integral 



Qvertex 



{2nf2^Mo 



A^^\p-q) 



{2ny 



/ 



d^k 



{{k + q + kiY — m?){{k + p — k^Y — m?)k'^ 



d'^kid'^k2d% 



g{ki)g{k2)g{k^). 



(440) 
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Of course, we choose the causal approach to the problem. The third order vertex contribution to the ^'-matrix 
can be written 

S^'^ = -e' J d^x.d^x^d^xs : cp''{x^)t^"'n^l, ^2, x^^x^) : A{xs)g{x^)g{x2)g{x^) + ... (441) 

The ^'-matrix element containing the external field is, correspondingly 

1 11 



Qvertex „3 



d'^xid'^X2d'^X3 / d'^kid'^k2d'^ks 



yielding 



(27r)32v^(27r)2(27r)6 
J d^k'Ts''''''''{xi,X2,X3)e^^'''e-''''''e-''''''''e-''''^^ 

1 



(442) 



-ivertex 



— — e 



(27r)327p^ 



d^kid%d%A'''\p -q- e(ki + k2 + ks)) 



g{eki)g{ek2)g{ek3), 



(443) 



where (note that tl^'^*^^{xi, 2:2, Xs) is translation invariant) 



ir''%P,q) = d'y,d%tr''%xr,X2,X3)e'Py^+"'y\ 



(444) 



Again, one can factor out the first order scattering matrix element, and using the abbreviations Pi — p — eki 
and qi — q + ek2 leads to 



Qvertex 0(1) 



d kid k2d k2.D{pi,qi)g{ekx)g{ek2)g{ek3). 



(445) 



In order to calculate t^^''*^^, we must first construct the causal distribution d\^'^*'^^{pi, qi), which is given by 

1 



d^k 



[D^n:\k - p,Wr:{k - q,)Dt\k) - Di;;\k - q,)bZ{k - pi)bt\k) 

+D^r:\k - q{)bz\k - p{)b^-\k) - b^r^\k - p,)bz\k - qi)bt\k) 

+b^^\k - qiM-\k - p,)b'Ak) - bit\k - p,)bl^ \k - q,)b^{k)]. 

The full calculation of d^^'^*^^ can be found in [7, 29]. The result is 



sgn(pio)©(p? - m^) log 



(446) 



4(27r)VA^ 

qf - - \ 



iV(l-^ 



Pi 



-sgn(gio)©(gi - m^) log 



p\ - m2 - - ^) - ViV(l - ^) 
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sgn(Po)0(^^ -4m2)log 



4m2 



p2 



(447) 



The time-ordered distribution is obtained from d"^^ by distribution splitting, i.e. from a subtracted 

dispersion integral according to the vertex scaling degree a; = 0. The infrared divergence is contained in the 
two first logarithmic terms of Eq. (447), and one may write 



J. 



div 



2^ 



f hiviti 



Idivit pi,tqi) 
t + iO) 



(448) 



with a first term 

Jdivi — 



8VA^(27r)6 



sgn{t)e{tV, - m') 



^^pm-m'' + f'{ql-p^q^) + f^ 



N 



(449) 

This integral can be evaluated in a straightforward manner and leads to expressions containing Spence func- 
tions and logarithms. We restrict ourselves to the term which contains the infrared divergence 



J. 



divi 



8\/iV(27r)6 
log 



log 



2 2 

m — PI 



pI 



log 



^pigi -m^ + ql- p^q, + ^{l - ^ 



^pigi -w? + ql - p,q, - ^/N{1 - ^) 



8VA^(27r)6 



pI 



log 



log 



Piqi - VN 



2 2 

m — PI 



pI 



log 
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Now we use the explicit form of pi and qi 
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The divergent vertex contribution to the cross section follows 
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Obviously, the self energy and vertex infrared divergences now cancel the divergence generated by the bremsstrahlung 
process: 
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This shows that the adiabatic limit g ^ I exists in the causal formalism for the inclusive cross-section. For a 
discussion of the uniqueness of the adiabatic limit we refer to [7, 30]. 

We conclude this section by highlighting the qualitative picture of the calculations given above. The switch- 
ing of the interaction with a test function g, which vanishes for large space and time distances, corresponds 
to a gedanken experiment where the charged mesons are liberated from their (scalar) "electromagnetic" field. 
The non-perturbative description of interacting fields is highly non-trivial and a hitherto unsolved problem. 
The good news is the fact that physical observables can be constructed in our model theory in an unambiguous 
way in the limit g ^ I, where the interaction becomes permanent. 



7 Gauge Theories 

7.1 Spin 1 

As we have emphasized before all knows interactions in nature can be described by quantum gauge theories. 
Gravity can be described within a very similar causal setting as "ordinary" spin-1 gauge theory, as will be 



58 



shown below. However, higher order perturbative quantum gravity holds the highly non-trivial problem of 
non-renormalizability, which may potentially show up in the causal framework as a violation of perturbative 
quantum gauge invariance. At least, the theory is still consistent at second order in the gravitational constant 
and may provide an effective description of the interaction. 

By quantum gauge theory we mean a theory which has a gauge invariant S-matrix. This is different 
from classical gauge invariance where the classical Lagrangian is gauge invariant. Instead we define gauge 
invariance for the time-ordered products constructed by causal perturbation theory as described in sect. 3 
using the gauge variations of the free fields. 

One is tempted to define perturbative gauge invariance simply by dgTn — 0, but this is not correct. To 
find the right definition let us consider QED where we certainly know what gauge invariance means. Ordinary 
spinor quantum electrodynamics is constructed from 

Ti(x) = ie : ^{x)-f^''ijj(x) : ^^(x). (461) 

The free Dirac fields ijj, ip have zero gauge variation dqip = = dqip, but dqA^ = id^u. Then we obtain 

dqTi — — e : ■07'*'0 : d^u — ie9^(i : ■07'*'0 : u). 

Here we have used current conservation 

which follows from the free Dirac equations. We see that dqTi is not zero, but a divergence 

dqT, = ia^Tf/,, (462) 

where 

Tf/^ = ie : ^-f^ip : u (463) 

is called Q-vertex in the following and Eq. (462) establishes the first order gauge invariance. 

It is not hard to generalize this to higher orders. If we freely interchange dq and the time ordering we can 
write 

dqTn = dqT{Ti{xi) ■ Ti(a;„)} 



= T{Ti{xi) . . . dqTiixi) Ti{xn)} 



1=1 

n 



= J2 T{Ti{xi) . . . td^T,),{xi) Ti{xn)}. (464) 
1=1 

The time ordered products herein have to be constructed correctly by the causal method, using the Q-vertex 
from Eq. (463) at xi instead of the ordinary QED vertex Eq. (461). Again formally taking the derivative out 
of the T-product we get 

n ^ 
O 



dqTn = i ^T{Ti{xi) . . . T'^/^{xi) ■ ... Ti{xn)} 
1=1 ^' 



1=1 <■ 

This equation certainly holds for xj ^ Xk, for all j ^ k, because there we can calculate with the T-product 
in the same way as with an ordinary product. But the extension to the diagonal xi = . . . = Xn produces local 
terms in general, both in Eq. (464) and in Eq. (465). If it is possible to absorb such local terms by suitable 
normalization of the distributions T„ and T^^^ then we call the theory gauge invariant to n-th order. We want 
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to emphasize that perturbative gauge invariance not only means that dqT^ is a divergence, the divergence must 
also be of the specific form Eq. (465) involving the Q-vertex. 

Now we check what perturbative gauge invariance defined by Eq. (465) means for the total S-matrix. 
Applying the gauge variation dq to the formal power series we obtain 

oo . » n 

dQS{g) -Y.h ^'^^ ■ ■ ■ E(^^'^n/J^(^l) ■ ■ ■ 9{Xn). 

n=l 1=1 

Since the test function g{x) is in Schwartz space, we can integrate by parts 

^~^Ia I ^^^^ ■ ■ ■ ^n/ldM ■ ■ ■ ■ ■ ■ 9{Xn)- 

n=l 1=1 

If it is possible to perform the so-called adiabatic limit g ^ 1 here, then the right-hand side goes to zero and 
we get the naive definition of gauge invariance of the S-matrix 

\imdQS{g)^Q. (466) 

The adiabatic limit exists if all gauge fields are massive. It does not exist for the time-ordered products if some 
gauge field is massless. In this case Eq. (466) is meaningless and we must use the perturbative definition in 
Eq. (465). The latter is really at the heart of gauge theory because it determines the possible couplings Ti. 
This we are now going to show for massless spin-1 gauge fields. 

We consider a collection of vector fields A^(x), a = 1, . . . and ghost fields Ua{x) with anti-ghost fields 
Ua{x) quantized in the usual manner according to Eq. (156). It is our goal to find all possible gauge invariant 
self-couplings Ti [x) of these fields. Since the gauge variation dq generates ghost fields from vector fields, 
it is pretty clear that Ti must involve ghost and anti-ghost fields as well. But we assume ghost number =0, 
so that u and u must appear in pairs. It is sufficient to consider triUnear couplings proportional to a coupling 
constant g, quadrilinear ones proportional to g'^ correspond to T2 and should come out automatically in the 
causal construction. We therefore start from the following general ansatz [31] 

Ti{x) = ig{fl,, : A,,A^,d^A>l : +/1 : A^^A^d^ A^, : 

+fL ■ A^aUbd^u, + /I : {d^^A.^aWi, : +/1 : A^a{d^u,)u, :}. (467) 

Here we have further assumed that Ti is a Lorentz scalar, and for the sake of simplicity we only consider CP 
conserving terms here. Ti being a Lorentz scalar, we need an odd number of derivatives in each term. We only 
consider one derivative because with three the theory is not renormalizable. The /^^^ are arbitrary constants, 
but unitarity requires a skew-adjoint Ti 

Tl{x) = -T,{x), 

so that the /'s and g must be real. This was the reason for the imaginary i in Eq. (467). Since the Wick 
monomial in the second term is symmetric in a and h, we assume 

J abc J bac 

without loss of generality. The reader easily convinces himself that there is no further possibility to contract 
the Lorentz indices and place the derivative. All double indices including a, 6, c are summed over. 
Next we calculate the gauge variation 

dqT^ = + A^^d^u.d-A'^^ + A^^^.d^d'^u,) 

+fl,{2d,UaAld,A^, + A,,Ald,d^u,) 
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+f!,,{d^Uad''ubUc + A^ad''u,d,A:)y (468) 

The last temi in the second and the first one in the fourth line vanish due to the wave equation. To simplify the 
notation we do not write the double dots for normal ordering anymore, all products of field operators with the 
same argument are normally ordered if nothing else is said. 

For gauge invariance the expression Eq. (468) must be a divergence 

We therefore write down a general ansatz for Tj^^^ as well: 

i^, = g[gl,^d''uaA,,A''^ + g^^UaA^.d^ A^ 

'^dlbc'^i^'^aA'j^A'^ + g'^f^^UaduA'^A'^ 

+gl,,UaA-^d,A^^ + gl,,u,u,d^u, + gl,p^u,u^u}^. (469) 
The symmetry in the first and antisymmetry in the sixth term give the relations 

9ahc 9achi 9abc 9bac' 

This ansatz for T^^-^ can be further restricted using the nilpotence property 

idQd^^, = d^Tr = 0. 

Substituting Eq. (469) and collecting terms with the same field operators we obtain the following homogeneous 
relations: 

d'^Uad^d^UbA"^ : 2glt,^ + g^f, + gl,^ - gl^ - gl, + ^^.^ = (470) 

d'^Uad.Ald^u, : 2gl,^ + gl, - g%^ - g^, = (471) 

d^Uad'^u.d^A'', : gl^ - gl, + g^^ - gl, + g^, - gl, = (472) 

Uad^u.d.d^A'^, : ^1 + 9lbe + albo - ale - ale = O (473) 

Uad>'d,u,d^Al : gl,^ + gl^, + gl,^ + gl^, = 0. (474) 
First order gauge invariance according to Eq. (468) now implies linear relations between the /'s and g's: 

d^u^d^A^.A'^^ : -fl, = 2gl,^ + gl, (475) 

d'^u,A,,d,A: : -2/2,^ - = + g^, (476) 

Uad^A^d.A^^ : -f^,, - /i, = gt,, + g^, {All) 

d.UaUbd^Uc-. -/l = 2^L + (/L (478) 

d.Uad^'UbUc : -f!bc + fbac = " dac (479) 

d.d^UaA'^A^, : -/I - = + 9lb (480) 

d.Uad.A'^A:: -fL = gL + 9Lb (481) 

Uad,A^,d,A^, : gl, = -gl, (482) 
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Uad.d^AlA':^ : -fl, = gl, + gl,^ (483) 

Uad^A^d^A^, : gl,, = -^L,. (484) 

All information comes out of this linear system. Since the elimination process is somewhat tedious, we 
give all details to save the readers time. Let us interchange h and c in Eq. (481) 

-/a\c = 9lc, + 9lc, (485) 

and add this to Eq. (481) 

-flbc - fLb = gibe + glcb + gibe + glcb- (486) 

By Eq. (484) g^ drops out and by Eq. (480) the right side is equal to 

= -fL - fL- (487) 

This implies 

fabc ~ fcba = fbca ~ facb- (488) 

Let us now decompose /^^^ into symmetric and antisymmetric parts in the first and third indices: 

fabc — ^abc ~l~ fabct dabc — dcbat fabc — fcbai (489) 

then Eq. (488) implies 

fabc facb fcba fcab fbca fbac- (490) 

So we arrive at the important result that fabc is totally antisymmetric. The Jacobi identity follows from second 
order gauge invariance, hence, fate can be regarded as structure constants of a real Lie algebra. 

The total antisymmetry of fabc implies the total symmetry of dabc- Next we use the representation Eq. (489) 
in Eq. (475): 

-fcab-dcab^2gl^ + gl,^. (491) 
Here g^^^ is antisymmetric in 6, c according to Eq. (482), so that g^^^^ must be symmetric, hence 

gibe = - h dcab = - h dabc (492) 

gibe = fcab = fabc- (493) 

Now we can write Eqns. (480) and (481) in the form 

gibe + glcb = -'^daa = -2dabc (494) 
-facb - dacb = gibe + glcb- (495) 

Since g^^^^ is antisymmetric in 6, c according to Eq. (484), the symmetric part of this equation agrees with Eq. 
(494) and the antisymmetric part is given by 

-facb = h (gibe - glcb) + glcb- (496) 

Hence, we find 

gibe = gibe - fabc + dabc: (497) 

where Eq. (494) has been taken into account. 

Now we turn to Eq. (473) and substitute g^ from Eq. (497) 

gtbc = glab - 2^ac6 " (glcb - facb + dabc) - (498) 
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Using this in Eq. (472) we see that and cancel out so that finally 

9l,c = l/abc (499) 

Then Eq. (498) can be simplified to 

albc = slab - alcb - dabc- (500) 

Substituting this into Eq. (476) gives 

fabc--h(fL + dac-dabc). (501) 

follows from Eq. (478): 

fabc — ~'^9abc ~ 9abc — ~9abc ~ fabc- (502) 

On the other hand, from Eq. (483) we get a different result 

fabc = -9lba + fabc: (503) 

which implies 

9lbc = 9lba - Vabc. (504) 

Finally, from Eq. (477) we conclude 

fbac + 9bac = ~ ftab ~ 9lab (505) 

and Eq. (479) gives another symmetry relation 

f!bc + albc = fL + glac- (506) 

It is easily checked that with the results just obtained all equations are identically satisfied. 
Summing up we have obtained the following form of the trilinear coupling 

T, = ig[{fabc + dabc)A,aA.,d''A'^^ - lifl^ + gl^ - da.M.aA^.d^A-'^ 

-{9lbc + fabc)A^aUbd^Uc + ft,c{d^ A^a)ubU, + fl,,A^a{d>'u^)uc] ■ (507) 
The terms proportional to dahc give a divergence 

dabc{A^.aA,,d^A''^ + lA^aA'^d^A'^) = Idabcd" (A^aA'^.A^c) ■ (508) 

This can be left out because it does not change the S-matrix. Next it is important to remember the relation Eq. 
(505) which shows the antisymmetry with respect to b and c. Therefore, we have 

Cl = i9{fabc + 9lbc){d^A^a)UbUc = -\9{ftbc + 9lbc)dQ{UaUbUc)- (509) 

Such a term which is dg of "something" is called a coboundary in cohomology theory [32]. Using Eq. (509) 
in Eq. (508) we have to add the term with which is taken into account as follows 

-gl^^{A^aUbd^Uc + d^A^aUbUc) 
= -9lbcd^iA^aUbUc) + glbcA„ad^UbUc). 

Now Ti assumes the following form 

Ti = ig\ fabcA^aAybd^A'^ — \ {flbc + 9lbc)AiJ,aA'^d''A'^ 



63 



+ ^-gdai>cd''{A^,aA''^A^^) - iggl^^df^iA^aUbUc) + Ci. (510) 
Due to Eq. (506) the second and fourth term together give a second coboundary 

C2 = ig{fL + glc){A^ad^ut,u, - lA^aA^.d^Al) 

-\9{flbc + 9lJdQ{A,,Alu,). (511) 

The coboundary terms lead to an equivalent S-matrix as well. 

Omitting the trivial divergence and coboundary terms we arrive at the following final result 

Ti = igfabM^aA^kd'^A''^ - A^aUbd'^Uc). (512) 

This is the well-known Yang-Mills plus ghost coupling to lowest order. At second order, gauge invariance 
gives the remaining coupling terms of pure Yang-Mills theory (see [6], sect. 3.4). 

The real strength of the method comes out in massive gauge theories. Since in S-matrix theory the asymp- 
totic free fields are the basic objects, one has to start with massive gauge fields from the beginning. Then gauge 
invariance of first and second order has to work and fixes all couplings. In particular, a physical scalar field, 
the Higgs field is necessary to satisfy second order gauge invariance. But the Brout-Englert-Higgs mechanism 
and spontaneous symmetry breaking plays no immediate role in such an approach. For details we refer to [6]. 

In order to motivate the formal accomplishments constructed so far, we conclude by giving a rather quali- 
tative comparison of the present formalism to the texbook literature. Above, we observed that QED is gauge 
invariant, but the true importance of gauge invariance is the fact that it allows to prove on a formal level the 
unitarity of the S-matrix on the physical subspace [41]. The presence of a skew-adjoint operator in the first 
order interaction or the presence of unphysical longitudinal and timelike photon states causes the S-matrix to 
be non-unitary on the full Fock space, but it is on the physical subspace. In QED, ghosts are introduced only 
as a formal tool, since they 'blow up' the Fock space unnecessarily, and they do not interact with the electrons 
and photons. But in QCD, the situation is not so trivial, due to the self-coupling of the gauge fields. 

The gluon vector potential can be represented by the traceless Hermitian 3x3 standard Gell-Mann matrices 
A«, a = 1,...8 

^M = E^My=^^;:y- (513) 

a=l 

The A's satisfy the commutation and normalization relations 



rA" A* 



c 



lU^j, tr(A"A^) = 25„;., (514) 



and the numerical values of the structure constants fate = —fbac = —fact can be found in numerous QCD 
textbooks. Since we are working with a fixed matrix representation, we do not care whether the color indices 
are upper or lower indices. 

The natural generalization of the QED Lagrangian to the Lagrangian of purely gluonic QCD is 

^gluon = --^^T^GiJ.uG'^" = -^G^^u^a"^ (515) 

with 

G^u = d^A, - d,A^ - ig[A^, A,] (516) 
or, using the first relation of Eq. (514) 

G% = d,Al - d^A; + gfabcA'^Al. (517) 
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It is an important detail that we are working with interacting classical fields here, therefore the correspond- 
ing field strength tensor G contains a term proportional to the coupling constant in contrast to the free fields 
F/;"^ = df.Afj^^ - d^Afj^^ used throughout this paper. Cgium is invariant under classical local gauge trans- 
formations 

A^{x) ^ U{x)A^{x)U-\x) + -U{x)d^U-\x), (518) 

where U{x) e 5[/(3). 

We extract now the first order gluon coupling from the Lagrangian. The Lagrangian 

C-giuon = -\[d,Al - d^Al + gUcA^AlWd^Al - A^^ + gj,,,,AlAl\, (519) 

contains obviously the free field part (this terminology is not really correct, since we are dealing with interact- 
ing fields here) 

4«on --\ [Q.K - {d^^Al - a^^^] (520) 
and the first order interaction part is given by 

l^tan = - d.A'^M'^^'o'K'K^ - \9fa,cA\A^AWK ' A^] (521) 

= -pabcA'^Atid'^A: - d-^A^] = - aM^f] (522) 

= gUcA^^Ald'^A^^. (523) 

The first interaction terms comes out from classical symmetry considerations here; in the framework presented 
in this paper, it is the consequence of purely quantum mechanical considerations. 

Since we are working in Feynman gauge, we add the corresponding gauge fixing term Cgf to the La- 
grangian. Additionally, we add the ghost term which describes the ghost interaction. The total Lagrangian 
then reads 

^QCD — ^gluon + ^gf + ^ ghost 

(524) 

= Cgiuon - ^{^^,A^f + d^'uidf.Ua - gfabcUbA^^). (525) 

The classical ghosts are anticommuting Grassmann numbers, i.e. v? = v? = 0, uu = —uu. 
The BRST transformation is defined by 

6AI = i\{d^Ua - gfabcUbA^c), (526) 
5ua = -iXd^A'^,, (527) 

5Ua = pfabcUbUc, (528) 

where A is a space-time independent anticommuting Grassmann variable. The special property of the BRST 
transformation is the fact that the actions 

= J d^x Sg, + Sgnost = / d'x {Cgf + Cg^ost) (529) 

and Stotai = Sgiuon + Sgf + Sghost are all invariant under the transformation: 

5 S gluon = 0, 6{Sgf + Sghost) = 0. (530) 

The similarity of free quantum gauge transformation introduced in this paper to the BRST transformation is 
obvious. One important difference is the absence of interaction terms ~ g. Furthermore, the free quantum 
gauge transformation is a transformation of free quantum fields, whereas the BRST transformation is a trans- 
formation of classical fields, which enter in path integrals when the theory is quantized. Finally, the free gauge 
transformation leaves the T„'s invariant up to divergences, whereas the BRST transformation is a symmetry 
of the full QCD Lagrangian. How the two symmetries are intertwined perturbatively is explained in [35]. A 
more rigorous axiomatic approach is discussed in [36, 37]. 
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7.2 Spin 2 

The crucial test of the gauge principle is spin 2 where it should lead to a quantum theory of gravity. In this 
case we supplement the gauge invariance condition 

[Q, T{x)] = dgTix) = id^T^ix), (531) 

in the following way. Since dq and the space-time derivative da commute it follows from nilpotency that 

dadqT- = 0. (532) 

If the appropriate form of the Poincare lemma is true, this implies 

dgT^^ [g,r"] = i9^r"^ (533) 

with antisymmetric T"'^. In the same way we get 

[Q,r"^] =ia^r"^T... (534) 

with totally antisymmetric T"^'"' and so on. These are the so-called descent equations (similar to Wess-Zumino 
consistency conditions). It is our aim to find a solution of these equations describing the self-coupling of the 
symmetric tensor field h^'^ considered in sect. 3.1, Eq. (176). We recall the gauge variations Eq. (180): 

dgh"" = -'-{d'u" + d^u" - v'^'dau") (535) 
dQv!" = 

dqu*" = (536) 

where we now denote the Minkowski tensor by rf^ to distinguish it from Einstein's g^^ . 

The descent procedure starts from T'^^'^ which must contain three ghost fields u and two derivatives and is 
totally antisymmetric. To exclude trivial couplings we require that it does not contain a co-boundary dqB for 
some B ^Q. Then there are the following two possibilities only: 

d^u'^u^'d^.u^, d'^u^'dyv?. (537) 

Therefore we start the descent procedure with the expression 

+d''u<u^d^u^ + ^^'u^u^'^^,u") + a2{d''u>'d^u^u' - d^u^d^u'^u' 

-d^u^'df.u^u'' - d^u^'d^ifu^ + d^'u^'d^u'^u^ + d'^u^'d^u^u''). (538) 

Next we have to compute d^T"^'^ and this is equal to —idqT'^f^ by Eq. (534). To determine T"'^ requires 
an "integration" dg^. As always in calculus this integration can be achieved by making a suitable ansatz for 
T"'^ and fixing the free parameters. The following 5 parameter expression will do: 

+^d^u''d^u^hf"' + hdXdi.Wh^'' - (a ^ p). (539) 
Substituting this into Eq. (534) leads to 

h\ — — 2ai, 62 = — 2a2 = — 2ai, 63 = 2ai, 64 = — 4ai, 65 = —2a\. 
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An overall factor is arbitrary, we take ai — —1 which gives 

T"^ = 2{ui'd^u^d^h'^'' + u'^d^u'^d^h^'' - u^'d^u^'d^h^'' + d^u^'d^u^h^" + d^u'^d^u^'h^^) -{a^ (5). (540) 

In a similar way we compute dpT'^^ and make an ansatz for T". The latter now has to contain ghost- 
antighost couplings also. The precise form can be taken from the following final result: 

-2v^d^u^d^fi!' + 2dXd^u''u^ + 2u^d^dXu'' - 2v^d^d^u^u\ (541) 
The last step calculating d^T^ and setting it equal to —idqT gives the trilinear coupling of massless gravity 

-4u^'dpu^df,h''f^ + Ad^uf^d^uph^"" - Ad^u^'d^u'^h^f' + Ad^u^d^upK"^ . (542) 

This coupling should have something to do with general relativity. To see this we leave quantum field 
theory aside and take the metric tensor g^^ as the fundamental classical field. The indices are no longer 
Lorentz indices, they are raised and lowered with g'^^ itself which is defined as the inverse g^i^g^'' — (^;^- One 
also introduces the determinant 

g^deig^,^. (543) 
Our starting point is the Einstein-Hilbert action given by 



Seh^-^J d^'x V^R, = 327iG, 



(544) 



where G is Newton's constant. R is the scalar curvature 

R = g^'^'R,.. (545) 

which follows from the Ricci tensor 

R/j-u — da^'^i, — dui^ila + r^/jF^^ — r°^rf^, (546) 

where 

= h9'"''{gp^^,'y +5m7>/3 -gpjw ) (547) 

are the Christoffel symbols. 

The variation of Eq. (544) is given by 

SehIq + ef] - SehIq] ^ejd'x [^^^gg-^^R,pr{x) 

+ j d''x^g''('(R^p[g + ef]-Ro^p[g\)+0{e'). (548) 
By calculating in geodesic coordinates one finds that the last term vanishes. Since 

d ^ 1 
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,al3 _ ^ al3 
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-V^{-l9,ug''^ + S';:5^), (549) 

we finally obtain 

SEH[9 + ef]-SEH[9]^e J d^x y/^g(-lg^^R + R^^^P^ix) + 0{e''). (550) 

This implies Einstein's field equations in vacuum 

Rap - = 0. (551) 

For this reason the Lagrangian 

Leh = -\V^R (552) 

can be taken as starting point of the classical theory. 

A glance at Eq. (546) and Eq. (547) shows that the first two terms in Eq. (546) contain second derivatives 
of the fundamental tensor field g^j^i,. This defect can be removed by splitting off a divergence. We rewrite the 
first term in Eq. (546) as 

v/=^^^T^,,« = (v^^'^'^r^J,, _r-^(^^M.),^ (553) 
and calculate the last derivative with the help of 

5'^^a=-^^a/'^-^^/J/^ (554) 

Proceeding with the second term in the same way we find 

^gR = - (v^^'^T^, - yf^gg^'^Tl^) (555) 

where 

G = g^\T':,T^,o.-njip)- (556) 
Since the divergence in Eq. (555) does not matter in the variational principle, we can go on with the Lagrangian 

2 



gg'^irZ^Tl^-TlJi^), (557) 



which contains first derivatives of g only. 

For the following it is convenient to remove the square root by introducing the so-called Goldberg 
variables 

r = V^99^\ ~9,u = {-9)-'''g,u. (558) 



Using 



dpg = 99a09"'^ 



9,9^'' = {-9)-'^'{r:p-hr''~9a(3r^p): 

dpg„u = {\~gnu~gapg'^'^,p -9noc~9upT'^ ,p ) (559) 



in Eq. (547) we obtain 



= \\\^T9^^r\',^\^~9^.r\p-~9p^r^.', 
-~9,,r^.p +r%,~9pur.p - \r'~9p^~9,.r\p ) • (560) 
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This enables us to express the Lagrangian L in Eq. (557) by g^i,. It is simple to compute the second term 

= -ir'',,~9.f9'^. . (561) 
But the first term in Eq. (557) requires the collection of many terms, until one arrives at the simple result 

-gapgpar^r^^. r'^.v +hr^~9^ur'',a -gpar^p ) ■ (562) 

Then the total Lagrangian is given by 



9i.ur\a9p.r'',/3r'')- (563) 



1 

~4 

To make contact with quantum field theory on Minkowski space we consider the situation in scattering 
theory where at large distances in space and time the geometry is flat and given by the Minkowski metric rj'^". 
Then we write the metric tensor as a sum 

g>'''{x) =7]'"' + Kh>"'{x). (564) 

We do not assume that the new dynamical field h^^{x) is small in some sense, it only goes to zero at large 
distances because of the asymptotically flat situation. The indices of h^^^ are ordinary Lorentz indices which 
can be raised and lowered with the Minkowski metric. Then the inverse of Eq. (564) is given by 

g^vix) = r]tj,,, - Khf,„{x) + K^hf.ah" ^- .... (565) 

Substituting these expressions into Eq. (563), the Lagrangian L becomes an infinite sum 



L = ^«;"L("). (566) 



n=0 



Here is the proliferation of couplings which is typical for gravity. It can be traced back to the infinite series 
Eq. (565). The three terms in Eq. (563) give the contributions 

Li = [-Vap + f^Kp - fi'KXp + • • • )^"^- h"'',^ (567) 

L2 = k (Vap - l^hap + n^haa'h'^' " • • • ) (^/3a " l^h/3a + /i/3/3' /if " • • • ) 

X {rj^" + Kh^'')hP^,^ h""",^ (568) 
X (r?"'^ + Kh'"^)hP'',a hP%^ . (569) 

The lowest order 

= ^h-^,^h':^ - h'^f'„h^^,f,-^h,^h-, (570) 
where h — h'^^j^, defines the free theory. Indeed, the corresponding Euler-Lagrange equations reads 

□ _ i^«/3g/^ _ /^^M,/?^^ -/l'^'^'",^ = 0. (571) 
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Both Eq. (570) and Eq. (571) are invariant under the classical gauge transformation 

The gauge can be specified by the Hilbert condition 

/i°'^,;3= 0. (573) 
This can be achieved by choosing the solution of the inhomogeneous wave equation 

as gauge function in Eq. (572). In the Hilbert gauge the equation of motion Eq. (571) gets simplified 

n/i°^^ - I'q'^'^nh = 0. (574) 

Taking the trace we conclude 

Uh = 0, Uh"^ = 0, (575) 

so that we precisely arrive at the free tensor field as it was assumed in the QFT. 
The first order coupling 0{k) in Eqns. (567)-(569) can easily be computed 

The first three terms herein agree precisely with the first three terms in Eq. (542). The last two terms and the 
forth and fifth terms in Eq. (542) are divergences. This is due to Lorentz contraction of the two derivatives. 
Indeed, if fi, /2, /s are massless fields satisfying the wave equation, then the following identity is true 

2dJid'^hh = 5"(5a/l/2/3 + hdj2h - hf2dah)- (577) 

Since divergence couplings do not change the physics, the coupling Eq. (542) derived from spin-2 quantum 
gauge theory agrees with general relativity in lowest order. It agrees at higher orders, too (see [6], sect. 5.7). 
The gauge principle even works in massive gravity [38]. The cohomological nature of gauge invariance was 
analyzed in [39]. 

The approach presented above is perturbative in nature and lives on a trivial background. Presently, no fully 
satisfactory theory of quantum gravity exists, and other ambitious approaches like, e.g., loop quantum gravity 
aim at a formally background independent description of quantum gravity, and they are expected to give rise to 
spacetime itself at distances which are large compared to the Planck length. How Einstein's classical geometric 
view on spacetime is related to such a theory is another story. Here, we content ourselves with the observation 
that we have found a gauge principle which uses the cohomological formulation of gauge invariance in Eq. 
(531) etc for the time-ordered products, having the character of a universal principle. Consequently it must be 
respected in any conventional regularization method. 

8 Conclusion 

Causality is a fundamental guiding element for the construction of perturbative quantum field theories. Us- 
ing causality in conjunction with a proper mathematical handling of distribution theory enables one to avoid 
ultraviolet divergences in perturbative quantum field theory from the start. Whereas standard methods like 
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dimensional regularization have calculational advantages compared to the causal method, the causal method 
provides a mathematically well-defined construction scheme of the perturbative ^'-matrix. 

In this review, a condensed introduction and overview of the causal approach to regularization theory has 
been given, which goes back to a classic paper by Henry Epstein and Vladimir Jurko Glaser [10]. The causal 
approach was taken up by Michael Diitsch and Giinter Scharf in 1985. During the last two decades, several 
important aspects of the theory have been worked out, which constitute the basis of this review. It should be 
mentioned that several topics which are not part of this work have been treated in the recent literature, like e.g. 
interacting fields [40], a complete discussion of perturbative QCD to all orders was worked out [41], and gauge 
theories like the full standard model (including phenomena like spontaneous symmetry breaking) were studied 
in [42, 43]. Theories in dimensions other than four were also considered [44] and specific analytic calculations 
of multi-loop diagrams were carried out [45, 29]. Supersymmetric theories [46] have been investigated, and 
the causal method was generalized to field theories on curved space-times [47] and studied in the framework 
of light cone quantum field theory [48]. 

As mentioned before, there are severe conceptual differences between the causal method and other regular- 
ization methods, which make it difficult to compare the different approaches in a reasonable way. Therefore, 
specific examples have been used in this work in order to demonstrate the differences and connections between 
the causal and dimensional regularization. On the one hand, dimensional regularization has many attractive 
features concerning the preservation of gauge invariance and in actual calculations due to its well established 
methods. On the other hand, the causal method is a strictly mathematical approach without any "intuitive" 
aspects like continuous spacetime dimensions. Furthermore, the formulation of quantum gauge invariance 
found during the study of gauge theories in the causal framework has a cogent structure when compared to 
the standard BRST approach [33, 34]. In this sense, the causal method constitutes an independent framework 
in its own right with many attractive features. Critical issues like, e.g. axial anomalies can be discussed in a 
unambiguous manner, and the strong mathematical background of the method permits to apply it to problems 
on curved spacetime and to quantum gravity, as we have illustrated in the last section. 

The mere observation that ultraviolet divergences can be avoided by a proper mathematical construction 
of Feynman diagrams certainly puts some arguments in the literature concerning the short-range behavior of 
quantum field theories in connection with ultraviolet divergences into perspective. Several approaches to QFT 
have been developed so far, and it is obvious that all considerations presented in this review are based on 
perturbation theory. Even if perturbation theory is well-defined order by order, it is far from being clear that 
the perturbation expansion can be summed up for physically relevant theories, even Borel summability is most 
probably not fulfilled due to Landau ghosts and renormalons. Despite these problems, perturbative QFT is 
a very successful and promising approach, since theoretical predictions of physical quantities made by using 
renormalized Feynman graph calculations match experimental results with a vertiginous high precision. In 
all these calculations, one should not think that it is impossible to avoid ill-defined integrals, as the causal 
approach proves. However, it should be mentioned that formal infinities are admissible if they are treated 
within a rigorous mathematical framework. Recent developments by Alain Connes and Dirk Kreimer based 
on Hopf algebras [49] have lead to some profound understanding how to "absorb" ultraviolet divergences in a 
consistent manner by a redefinition of the parameters defining the QFT. The Hopf algebra approach has also 
been applied to the causal Epstein-Glaser approach in [50], in order to overcome the separation between the 
causal method and mainstream QFT. 

Richard Feynman in his Nobel lecture remarked: "I think that the renormalization theory is simply a way 
to sweep the difficulties of the divergencies of electrodynamics under the rug." This problem has been solved 
by the causal method, at least on a perturbative level. 
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A Special Distributions in 3+1-Dimensional Spacetime 

In this appendix, we give a condensed account of the most important properties of the causal commutators and 
propagators used in the present review. The distributions used in the causal framework typically differ from 
the most common definitions found in the literature and by a sign or a normalization factor, since we use the 
"mathematical" symmetric definition of the (inverse) Fourier transform. Accordingly, the distributions used in 
the present text are related to the distributions below by the simple redefinitions 

{'^^fD^,rem = -A^-*(A;), {27:fD'^^\k) = -A±(A;) (578) 

in momentum space and by Dp\^^{x) = —A^'"\x), Dm\x) = —A^{x) in real space, omitting potential 
mass indices. 

The free (non-interacting) neutral scalar quantum field (p{x) for particles with a given mass m is given by 

(kx = k^x^', k^ = y(P + m^)) 

1 r d^k 

^{x) = ^-{x) + ^\x) = j + a\k)e-^^'% (579) 

where ^p~{x) and (p^{x) refer to the corresponding frequency parts, respectively, whereas the charged field 
i^dx) has the form 

Mx) = ^:{x) + ^ti^) = J ^[«(^)^"''' + b\k)e+''% (580) 

The commutators of the operator- valued distributions a{k), b{k) ("annihilation operators") and a^(^), b^k) 
("creation operators") are 

[a{k), a^{k')] = [b{k), b^{k')] = S^^\k - k'), (581) 
[a{k),a{k')] = [b{k),b{k')] = [a\k),a\k')] = [b\k),b\k')] = V^,^', (582) 
and the annihilation operators destroy the unique perturbative vacuum |0) according toa(A;)|0) = b{k)\0) = 

v£ 

The commutation relations of the scalar fields lead to the so-called positive and negative frequency Jordan- 
Pauli distributions 

Aiix) = -tiv^ix), if^m = -tiOliif^ix), if^iom) , (583) 
with the distributional Fourier transforms 

A^(fc) = J d^x A^{x)e"''' = T{27Ti) e(±fc°)5(fc' - m'). (584) 

S is the one-dimensional Dirac distribution depending on k'^ — k^k^ — {k^Y — (k^Y — {k'^Y — {k^Y — 
The fact that the commutator 

[^{x), 99(0)] = lAiix) + iA;^{x) =: iAm{x) (585) 

vanishes for spacelike arguments < due to the requirement of microcausality, leads to the important 
property that the Jordan-Pauli distribution has causal support, i.e. it vanishes outside the closed forward 
and backward light-cone such that 

supp Am(x) CV'UV^ , ^ {x\x^ >0, ±x° > 0} (586) 

in the sense of distributions. 
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The retarded propagator A'!^^{x) is defined in configuration space by 

A-*(x) = e(x°)A^(x), (587) 
leading to the Fourier transformed expression 

(^) = T-o 1 , ./nn - (588) 

The Fejaiman propagator is given in configuration space by the vacuum expectation value 

A^{x) = -t{0\T{^,{x)^l{0)\0) = -t{0\T{^l{x)^M\Q) - t{0\T{^{xMO)\0), (589) 
the well-known distributional Fourier transform reads 

In the massless case, one has 

° " J (27r)4A;2 + ^0 " 4^^^^^ " A^^ x^ ~ 4^ ^ ^ 

where T is the time-ordering operator, P denotes principal value regularization. 

It is straightforward to show that the distributions introduced above fulfill the distributional differential 
equations displayed below. From the wave equation (□ + m'^)ip'^^\x) = (□ + rn'^)ip^{x) — follows 

(□ + m2)A^(x) = (□ + m2)A^(x) = 0. (592) 

Furthermore, one has 

(□ + m^)A^(x) = {d^d^ + m^)AlXx) = -5W(x), (593) 

and 

(□ + m^)A'-^{x) = -5(^)(x). (594) 
The Feynman propagator and the retarded propagator are related via 

^ret ^ + A" . (595) 
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